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ABSTRACT 


The  ROSE  balloon/AN-FPS-16  radar  system  is  designed  to  report 
small-scale  winds  and  shears  by  polynomial  smoothing  of  the  radar  tracking 
data.  The  most  successful  mathematical  technique  for  fitting  the  data  and 
computing  winds  and  shears  is  presented  along  with  other  techniques  which 
were  considered  but  ultimately  rejected. 

Analysis  of  the  self-induced  balloon  motion  is  attempted  by  classical 
physical  methods  and  by  a  Power  Spectral  Density  Technique. 

A  criterion  for  evaluating  the  effect  of  modification  to  the  standard  ROSE 
is  developed  and  applied  to  modified  ROSEs  and  the  newer  small  lightweight 
balloon. 
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INTRODUCTION 


The  standard  ROSE  balloon  is  a  metalized  nonexpanding  sphere  having  a 
diameter  of  two  meters.  When  tracked  by  AN/FPS-16  radar,  the  space-time 
coordinates  of  this  rising  balloon  can  be  determined  accurately  at  a 
frequency  of  10/sec. 

This  report  is  a  chronological  presentation  of  part  of  the  work  per¬ 
formed  under  AF  19(604)-7450  sponsored  by  AFCRL.  The  first  section 
describes  only  the  mathematics  which  were  developed  and  used  in  an  attempt 
to  accurately  describe  the  balloon's  velocity  and  acceleration.  As  the  reader 
will  note,  the  author  became  increasingly  aware  that  the  velocities  and 
accelerations  of  the  standard  ROSE  balloon  did  not  represent  the  wind.  In 
particular,  the  acceleration  of  the  balloon  in  no  way  represents  the  wind 
shear.  Consequently,  Sections  2  and  3  deal  with  the  unsuccessful  effort  to 
eliminate  "noise"  (regardless  of  cause)  from  the  balloon  motion  in  order  to 
find  the  proper  wind. 

Concurrent  with  the  work  described  in  Sections  1-3,  Murrow  and  Henry 
(working  independently)  determined  by  physical  experiments  that  the  standard 
ROSE  balloon  was  definitely  not  a  proper  wind  sensor--at  least  over  the 
first  200  feet  of  its  ascent.  Reid  of  AFCRL,  who  was  also  aware  of  this, 
showed  that  unpredictable  aerodynamic  forces  were  acting  on  the  standard 
ROSE  over  at  least  the  first  20,  000  feet  of  its  ascent.  Our  study,  especially 
the  work  described  in  Sections  2  and  3,  indicated  that  there  was  no  good  way 
to  mathematically  separate  the  aerodynamic  forces  acting  on  the  balloon 
from  those  induced  on  the  balloon  by  the  wind.  It  thus  became  apparent  to 
all  concerned  that  the  standard  ROSE  was  not  an  adequate  wind  sensor. 

A  physical  modification  of  the  standard  ROSE  balloon  was  therefore 
undertaken  by  AFCRL.  Such  work  had  already  been  started  by  Scoggins  of 
NASA.  His  modified  ROSEs  are  called  "  Jimspheres.  "  Reid  and  Lenhard  of 
AFCRL  devised  several  different  modifications  and  developed  some  smaller 
smooth  ROSEs.  The  description  of  these  balloons,  the  mathematical 
treatment  used  in  describing  their  motion,  and  the  statistics  employed  in 
comparing  the  modified  balloons  to  the  standard  and  to  the  smaller  ROSEs  is 
described  in  Sections  4,  5,  and  6. 

The  conclusion  reached  in  these  sections  is  that  the  modification  has 
improved  the  wind  sensitivity  and  will  permit  use  of  the  techniques  tor  pre¬ 
dicting  wind  and  wind  shear  which  are  described  in  Section  1. 


1. 


MATHEMATICAL  TECHNIQUES  FOR  COMPUTING  WINDS  AND  SHEARS 


That  the  ROSE  balloon  constituted  a  good  wind  sensing  system  was  the 
basic  assumption  for  the  various  mathematical  techniques  developed  for  com¬ 
puting  winds  and  shears.  Sections  1.1  and  1.2  describe  two  of  these  techniques 
and  the  empirical  basis  on  which  they  were  rejected.  Sections  1.3  and  1.4 
describe  the  other  techniques  which  seem  preferable.  Although  none  of  these 
methods  were  acceptable  prior  to  balloon  modification,  it  now  appears  that 
the  techniques  described  in  Sections  1.3  and  1.4  will  provide  valid  wind  and 
shear  values  when  applied  to  radar  tracking  data  from  the  modified  balloon. 

The  mathematical  rationale  upon  which  all  these  techniques  are  based 
is  as  follows.  The  standard  equations  of  motion  for  a  balloon  moving  in  the 
atmosphere  without  external  forces  are: 

Mx  =  1/2  p  Afi  Cd(v|  (x  -  Wx)  (1) 

My  =  1/2  p  Afi  CDlv!  (y  -  Wy)  (2) 

M'z  =  1/2  p  Afi  Cd  |  v|  (z  -  Wx)  +  (m  -  p  Vfi)  g,  (3) 


where  m  is  the  mass  of  the  balloon  and  M  is  the  total  effective  mass  (see 
Appendix  1).  Combining  the  first  and  third  equations  and  then  the  second  and 
third  gives  the  equations  for  the  component  wind  in  the  x  and  y  directions: 


W 

x 


X 


Mx  (y  -  Wy) 

Mj  ■  (m  ■  pVB|  g 


and 


W 

y 


My  (i  -  Wz) 
y  ’  Mz  -  (m  -  pVfi)  g 


(4a) 

(4b) 


The  second  term  in  both  equations  may  be  neglected  because  of  its 
insignificance  when  compared  with  the  first  term.  The  wind  equations 
modified  in  this  way  would  be 


W  =  x  and 
x 


(5a) 

(5b) 


Winds  computed  by  this  equation  are  known  as  point  winds.  The  average  wind 
in  a  layer  is  given  by  the  equations 


W 


x 


and 


X2  "  yl 


(6a) 

(6b) 


where  x- ,  y^,  and  L  are  the  space-time  coordinates. 

The  wind  shear  in  the  component  directions  is  defined  by  the  equations 


d<Wx> 

dh 


and 


d(W  ) 

y 

dh  * 


(7a) 

(7b) 


The  chain  rule  of  differentiation  gives 
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dt_ 

dh 


and 


dt_ 

dh 


(8a) 

(8b) 
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dt 


S 
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dh 


If  the  approximations  x 


and  z  s  h  are  made  the  equations  are: 


(9a) 

(9b) 


Quite  often  it  is  difficult  to  compute  the  second  derivative  accurately.  In 
such  cases  d(Wx)  is  replaced  by  AWX,  and  dh  is  replaced  by  Ah  which  are 
finite  differences.  The  component  snear  equations  then  become: 


(10a) 

(10b) 


1.  1  NINTH-DEGREE  POLYNOMIAL  METHOD 


Sx  = 


W  -  W 
x2  X1 

"V^I . 


W 


w 


and 


sy  = 


It  has  been  suggested  by  Scoggins  (May,  1963)  that  a  polynomial  approx¬ 
imation  to  the  FPS-16  radar  tracking  data  will  yield  valid  accelerations  and 
velocities  whirh  will  in  turn  give  realistic  values  of  winds  and  shears  from  a 
rising  sphere1.  It  is  most  economical  to  use  the  lowest  possible  degree 
polynomial  in  approximating  the  radar  tracking  data.  In  order  to  determine 
the  lowest  degree  necessary  the  data  was  approximated  by  a  zero  degree 
polynomial,  and  the  residual  was  compared  with  the  theoretical  tracking  error 
inherent  in  the  radar.  The  degree  of  the  polynomial  was  then  increased  until 
the  residual  was  less  than  the  tracking  error.  The  maximum  degree  used  was 
the  ninth  degree  even  though  the  residual  of  the  ninth  was  occasionally  larger 
than  the  tracking  error.  Since  the  radar  tracking  data  is  equally  spaced  at 
one-tenth  second,  orthogonal  polynomials  are  used  in  the  approximations.  The 
orthogonal  polynomials  are  especially  desirable  in  this  application  because 
they  allow  the  next  higher  power  term  to  be  added  without  recalculating  all 
the  previous  coefficients  as  is  the  case  with  ordinary  least-squares  poly¬ 
nomials.  Later  work  indicated  that  it  was  also  desirable  to  compute  the  anal¬ 
ysis  of  variance  (ANOVA)  table  to  determine  whether  the  higher  degree 
polynomial  terms  were  statistically  significant. 

1.1.1  ORTHOGONAL  POLYNOMIALS 


The  general  form  of  the  orthogonal  polynomial  for  n+1  data  points,  eval¬ 
uated  at  the  i  th  point,  is: 


=  A0P0,n(l)  +  AlPl,n(l)  +  +  AmPm,n*l) 


(H) 


where  i  takes  on  all  integral  values  from  0  to  n,  and  q  is  x,  y,  or  z  depending 
on  which  axis  is  being  considered  ind  m  is  the  degree  of  the  polynomial. 


A  is  the  coefficient  of  the  form 
m 


3 


A 


E  q  P  n(i) 

i=Q  1  rn>  n 
m  n  -> 

L  Pm 

uo  m'n 


where  is  the  i  £h  data  point. 


(12) 


^ ( i )  is  a  Legendre  polynomial  of  degree  m,  of  the  form 


m,  n 


<i>  =  £  (-l)k  (™)  <mkk>  “ 

k=0  k 


.W 


n 


Ck] 


where 


i  ^  =  i(i-l)(i-2)(i-3)«  •  *(i-k+ 1 ),  i^=  1  and  (k^  = 


m'. 


(m-k)ikl 


(13) 

(14a) 

(14b) 


The  evaluation  of  these  expressions  for  each  degree  from  zero  to  nine  is 
given  in  Appendix  6. 

The  Legendre  polynomials  were  generated  in  a  computer  program  instead 
of  being  obtained  from  standard  tables  so  that  transcription  errors  and  costly 
key  punching  might  be  avoided.  The  generating  program  was  written  in 
BALGOL  and  run  on  the  University  of  Dayton's  Burroughs  220  computer. 

Every  foreseeable  error-producing  operation  was  eliminated.  When  possible, 
integer  arithmetic  was  used.  When  it  was  not  possible  to  use  integer 
arithmetic,  18-digit,  double -precision  a  rithmetic  was  used.  The  orthogo¬ 
nality  of  each  set  of  numerical  values  of  Legendre  polynomials  was  checked 
using  this  property  of  orthogonality: 


n 

L 

i=0 


P  (i) 

m,  n 


0 


and 


n 

L 

i  =  0 


P  (i)  = 

m,  n 


m=2,  5 


m  =  7,  1  0 


(15) 


(16) 


The  values  of  the  coefficients  used  in  the  standard  reduction  program 
were  limited  to  8  digits;  therefore  the  results  of  the  generating  program  were 
truncated  after  the  8th  digit.  The  maximum  deviation  from  zero  occurring 
in  the  largest  sets  was  less  than  9  x  10  .  It  was  felt  that  these  coefficients 

would  produce  satisfactory  results  in  the  data  reduction  program. 


The  sums  of  squares  of  the  coefficients 


n 

L 

i  =  0 


we  re  also  com¬ 


puted  in  double  precision  and  included  at  the  end  of  the  list  of  each  set. 


A  sample  printout  of  one  set  of  coefficients  is  presented  in  Figure  1. 
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Figure  1.  Tabular  Listing  of  Orthogonal  Polynomial  Coefficients  for  25  Points 


The  physical  size  of  the  computer  used  in  the  computation  of  winds  and 
shears  imposes  a  practical  limit  on  the  maximum  number  of  Legendre 
polynomials  that  can  be  used  in  the  approximation  of  the  data.  The  maximum 
number  of  data  points  to  be  approximated  can  be  determined  by  considering 
the  vertical  velocity  of  the  balloon.  The  average  vertical  velocity  of  the 
balloon  is  about  20  feet  per  second,  producing  an  average  of  50  tenth-second 
data  points  per  100-foot  layer.  In  order  to  include  all  reasonable  vertical 
velocities,  the  lower  limit  was  set  at  10  feet  per  second  and  the  upper  limit 
at  40  feet  per  second.  This  corresponds  to  100  points  and  25  points  for  the 
upper  and  lower  limits  respectively.  The  Legendre  polynomials  were  com¬ 
puted  to  cover  the  range  of  25  to  100  data  points  per  fit  and  the  degrees  from 
one  to  ten. 

It  was  desirable  to  associate  an  altitude  in  the  layer  with  the  meteoro¬ 
logical  parameters.  The  altitude  chosen  was  the  mid-point  of  the  100-foot 
layer.  The  accelerations  and  velocities  were  evaluated  at  mid-point  of  the 
data,  and  these  values  were  associated  with  the  altitude  50  feet  above  the 
bottom  of  the  layer.  It  is  much  more  convenient  to  evaluate  the  polynomials 
at  one  of  their  original  data  than  between  two  of  them.  The  mid-point 
coincides  with  an  original  data  point  only  when  there  is  an  odd  number  of 
data  points  in  the  layer.  The  program  is  designed  to  assemble  an  odd  num¬ 
ber  of  data  points  over  the  nominal  100-foct  layer  which  might  be  as  much  as 
103  feet  in  depth.  Each  layer  was  considered  individually.  That  is,  there 
was  no  overlapping  of  data  from  layer  to  layer. 

1.1.2  TRACKING  ERRORS 

The  theoretical  tracking  errors  to  be  used  in  the  determination  of  the 
degree  of  polynomial  approximating  the  data  is  a  function  of  the  errors  in  the 
parameters  reported  by  the  radar,  namely  elevation  and  azimuth  angles  and 
the  slant  range  distance.  These  radar  parameter  errors  are  related  to  the 
space-time  coordinate  errors  by  the  geometry  of  the  coordinate  systems. 

The  space-time  coordinate  system  used  is  earth  fixed,  Euclidian 
3-space,  such  that 

i)  The  x-axis  is  oriented  north  -  south,  with  north  positive, 
li)  The  y-axis  is  oriented  east  -  west,  with  east  positive, 
in)  The  z-axis  is  perpendicular  to  the  x  -y  plane  with  upwards  positive, 
iv)  The  x  -  y  plane  is  tangent  to  the  earth  at  the  launch  site. 

The  equations  of  transformation  from  radar  parameters--Range  (R), 
Elevation  (E),  and  Azimuth  (A,  positive  clockwise  from  north)--to  a  Euclidian 
3  -  space  a  re: 


R 

cos 

E 

cos  A 

(17) 

R 

cos 

E 

sin  A 

(18) 

R 

s  in 

E 

(19) 
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Since  there  exists  an  error  inherent  in  the  radar,  the  errors  of  x,  y, 
and  z  with  respect  to  the  error  in  A,  E,  and  R  are  given  by 


D2  T  2.2 

*  R  L'e  ,,n 

d2  f  2.2 

=  R  sin 


_  2  ,  2 
E  cos  A  ♦  o  cos 
A 


2  Esin2A  j  +  aR 


2cos2£cos2A  (20) 


2  .  >  .  2  .  2  2  „  2. 
<r_  sin  Esin  A  +  <r  cos  Ecos  A 
E  A 


s2Esin2A  (21) 


2  2  2  _  t  2  .2  c, 

a?  s  R  o-£  cos  E  +  sin  E  I 

2  2  2 

where  a  ,  r  and  aD  are  the  errors  in  Azimuth,  Elevation,  and  Range 
A  £  K 

respectively  and  are  constant  with  respect  to  the  radar  system  tracking  the 
balloon. 


Since  the  data  is  given  in  x,  y,  and  z,  the  transformation  equation  can  be 

used  to  give  the  errors  <r  2f  2,  and  <r  2  in  terms  of  x,  y,  and  z  and  the 

x  y  « 

constant  radar  errors.  The  tracking  errors  in  this  more  useful  form  are: 


The  errors  must  be  computed  for  each  layer  since  the  erroro  are  functions  of 
the  x,  y,  and  z  values. 

1.1,3  DETERMINATION  OF  DEGREE  OF  THE  POLYNOMIAL 

As  outlined  in  a  previous  section,  the  lowest  degree  polynomial  possible 
in  approximating  the  radar  tracking  data  is  determined  by  comparing  the 
residual  of  the  m  th  degree  polynomial  fit  with  the  independent,  theoretical 
tracking  error  described  in  Section  1.1.2.  The  residual  of  the  zero-degree 
fit  is  compared  to  the  tracking  error;  if  the  tracking  error  is  larger,  a ’’good" 
fit  has  been  obtained,  in  this  case  the  mean.  If  the  tracking  error  is  smaller 
than  the  residual,  the  sum  of  squares  accounted  for  by  the  first-degree  fit  is 
subtracted;  this  new  residual  is  then  compared  with  the  tracking  error.  If 
the  residual  is  now  smaller  than  the  tracking  error,  the  linear  fit  is  "good"; 
if  the  residual  is  still  larger,  the  next  degree  (the  second)  must  be  added. 
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This  process  continues  until  the  residual  is  finally  smaller  than  the  tracking 
error,  or  until  the  ninth  degree  has  been  reached.  If  the  residual  is  still 
larger  than  the  theoretical  tracking  error,  the  process  is  discontinued,  and 
the  ninth-degree  polynomial  fit  is  accepted  as  being  "good." 


The  equation  of  the  residual  of  the  zero-degree  fit  is 


n 

£ 

i=0 


(qj  -  m2 
TH — 


(26) 


where  are  the  data  points  and  =  Aq  since  Py  n  (i)  =  1*  Rq  ^  is 
compared  to  cr^  (in  calculating  a-^,  the  mean  of  in  the  100  foot  layer  is 
used).  If  R0  is  less  than  or  equal  to  <r  a  zero-degree  polynomial  is  a 
"good"  fit.  If  Rn  is  greater  than  a  ,  a  higher  degree  fit  is  necessary. 

"» q  q 

The  expression  for  the  residual  of  the  j  +  1  st  degree  fit  is  given  by  the 
recurrence  relation 


R.  I  =  R.  -JiiUL 

j-n»q  j.q  n  +  i 


(27) 


where 


j+i.q2  =  l?o[Aj+lPj+1,n<‘)] 


=  Aj+i 


n 


i  2 


L  P 

i=0 


j+l,n 


(i) 


(28) 


The  coefficients  Aj,  A £,  ...»  Aj+j  are  of  course  available  from  the 
polynomial  without  further  computation. 

An  alternate  method  of  determining  the  lowest  degree  of  the  polynomial 
is  the  standard  statistical  technique  of  the  analysis-of-variance  table  (Crow, 
Davis,  Maxfield)  .  This  table  is  presented  as  Figure  2.  The  analysis  of 
variance  is  performed  to  test  the  significance  of  the  terms  AmPn^nU)  of  the 
polynomial.  First,  the  total  sum  of  squares  (which  is  the  sum  of  squared 
deviations  of  q  from  the  average  q)  is  computed.  As  shown  earlier  in  this 
section,  the  average  q  is  just  Aq,  the  first  term  of  the  polynomial.  Computed 
next  is  that  part  of  the  *otal  sum  of  squares  accounted  for  by  the  linear, 
quadratic,  cubic,  quartic  terms,  etc,  up  to  the  m  th  degree.  These  sums  of 
s qua  res, 


their  degrees  of  freedom,  and  the  corresponding  mean  squares  are  shown  in 
F  igure  2.  That  part  of  the  total  not  accounted  for  by  the  sums  of  squares  of 
the  various  degree  terms  is  the  residual  sum  of  squares.  An  "F"  value  is 
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computed  for  each  term  by  computing  the  ratio  of  the  mean  square  of  each 
term  to  the  mean  square  of  the  residual.  Each  computed  F  value  is  then 
compared  with  the  tabular  value  of  the  "F  Distribution"  at  the  5  percent  sig¬ 
nificance  level.  If  the  F  of  the  j  th  term  is  larger  than  the  corresponding 
tabular  value,  then  the  data  are  significantly  better  represented  by  including 
this  term  than  by  omitting  it. 

The  analysis -of- variance  criterion  was  not  applied  in  this  problem 
because  it  was  felt  that  the  independent  estimate  of  error  was  more  meaning¬ 
ful.  In  order  to  become  familiar  with  the  contribution  to  the  fit  of  each  term 
in  the  polynomial,  a  modified  form  of  this  table  was  computed  and  printed 
out  for  several  early  flights.  The  modification  was  to  reduce  the  size  of  the 
tables  by  printing  out  only  the  sum  of  squares  of  the  residual  after  each  term 
had  been  added  to  the  polynomial.  Some  "F"  values  were  then  computed  by 
hand  and  checked  for  significance. 


1.  1.4  COMPUTATION  OF  WINDS  AND  SHEARS 

1 .  1.4.  1  Winds 

There  are  two  ways  to  describe  the  wind  in  the  100-foot  layer;  1)  the 
average  wind  in  the  layer;  2)  the  wind  evaluated  at  a  given  altitude  known  as 
the  Point  Wind. 

1. 1.4.2  Average  Wind 

The  average  wind  in  the  layer  is  determined  by  a  simple  finite  difference 
between  the  highest  fitted  coordinate  point  in  the  layer  and  the  lowest  point 
divided  by  the  time  needed  to  traverse  the  layer.  In  terms  of  the  evaluated 
polynomial  fit,  the  average  wind  is  given  by: 


1.1.4.  3  Point  Wind 

Using  .the  equation  developed  in  Section  l,  the  wind  is  given  by 
WG  =  10  •  the  factor  of  10  is  necesrary  because  the  points  are  0.  1 

Mi  1  ^ 

second  apart  instead  of  1 . 0  second  apart.  Note  thatq^  is  the  first  derivative 
of  the  polynomial  evaluated  at  the  point  i.  The  mid-point  of  the  layer  was 
chosen  as  the  altitude  to  be  assigned  to  the  point  wind.  This  corresponds  to 
evaluating  the  polynomial  at  the  point  (i-n/2)  where  n  is  always  even. 

The  equation  for  the  first  derivative  is 


A,P'  Jr)  +  A\P\  (t)  +  ...  +  A.P'.  (?) 

1  l,nY'  3  3,n'2'  j  j,  m  2' 


(30) 


where  j  is  an  odd  number  either  equal  to  m  (if  m  is  odd)  or  m-1  (if  m  is  even). 
The  point  wind  is  then  given  by 


10 


=  10  •  3 


(31) 


W 

qn/2 


n  /2 


Obviously,  if  m  =  1  (the  polynomial  is  of  the  first  degree),  then  the 
point  \und,  W  ,  will  be  equal  to  the  average  wind,  W  . 

qn/2  ^  q 

The  derivatives  of  the  Legendre  polynomials  are  developed  and 
presented  in  Appendix  6. 

1 . 1.4.4  Shea  rs 


During  the  examination  of  the  results  of  preliminary  computations 
(particularly  the  ANOVA  tables)  it  was  decided  that  accurate  second  deriva¬ 
tives  of  the  polynomial  would  be  difficult  to  obtain,  especially  in  the  higher 
degrees.  Therefore,  the  simplified  form  of  the  shear  equations  was  chosen: 


S 

x 


W 

x 

IT 


S 

y 


(10a) 

(10b) 


It  was  thought  that  the  point  wind  (Section  1.1.4.  3)  was  a  better  representation 
of  the  wind  at  the  center  of  the  layer;  therefore  point  winds  were  used  in  the 
computation  of  the  shears.  The  altitude  layer  was  always  100  feet;  conse¬ 
quently  the  denominator,  h^  *  hj,  always  equals  100  feet. 

The  equations  for  shear  then  reduce  to 


TOT 


S  = 

y 


t  +  i 


■(yn/2) 


100 


(32a) 

(32b) 


where  the  subscript,  t,  denotes  any  given  layer  and  f  +  1,  the  layer  immediately 
above  it.  The  units  of  shear  are  sec"  . 

1.1.5  RESULTS  AND  CONCLUSIONS 

The  computer  program  was  written  to  compute  the  winds  and  shears  as 
outlined  in  the  previous  sections.  The  program  was  written  in  FORTRAN  II 
for  the  IBM  7094  computer,  checked  out  at  the  W right-Patte rson  A.F.B. 
facility,  then  forwarded  to  Eglin  A.F.B.  for  operation  on  their  system. 

Through  the  very  fine  cooperation  of  PGVMS,  Eglin  A.  F.  B.  ,  Florida  the 
program  with  all  of  the  polynomial  coefficients  was  implemented.  A  test 
ROSE  balloon  was  launched  and  tracked,  the  data  were  reduced  in  manner 
that  has  been  outlined,  and  the  results  were  forwarded  to  the  University  of 
Dayton. 

The  results  of  this  test  run  were  quite  surprising  in  that  while  the 
average  winds  looked  reasonable,  the  point  winds  (and  therefore  the  shears) 
were  quite  erratic  and  did  not  agree  with  the  average  winds.  A  sample  of 
this  output  is  shown  in  Figure  3.  As  indicated  on  the  tabulation  shown  in 
Figure  3,  the  output  parameters  are: 


1 1 


AIT!  luOt  »INC  CAST  NORTH  TOTAL  ANG  *  PT  SMfcAR  IASI  NORTH  TOTAL 

IOSO.  0.00  -9.03  5.3S  C  S 
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Figure  3.  Sample  Printout  of  Ninth-Degree  Method 


Altitude 
Wind  East 
Wind  North 
Wind  Total 
Angle 
M 

Shear  East 
Shear  North 
Shear  Total 
Point  Wind  East 
Point  Wind  North 


ft 

ft/sec 

ft/sec 

Knots 

degrees 

Highest  degree  polynomial  used 

Sec"1 

Sec"  1 

Sec"1 

ft/ sec 

ft/sec 


As  will  be  noted,  the  highest  degree  used  by  either  component  is  quite 
often  the  ninth,  which  was  the  maximum.  The  component  data  were  plotted 
vs.  time  and  although  ragged,  they  resembled  a  parabola--indicating  a  second- 
degree  polynomial  fit.  It  was  then  decided  to  add  monitor  variables  to  the 
printout  to  check  intermediate  computations  and  to  add  the  Residuals  of  the 
fits  as  successively  higher  degree  terms  were  added.  A  sample  of  another 
test  flight  is  presented  in  Figure  4.  The  monitor  parameters,  indicated  by 
the  +,  are 


ft  (average  space  coordinate) 
Tracking  error  for  this  layer 
Residual  for  zero  degree 
Residual  for  first  degree 
Residual  for  second  degree 
Residual  for  third  degree 
Residual  for  fifth  degree 
Residual  for  seventh  degree 
Residual  for  ninth  degree 


The  monitor  parameters  are  always  arranged:  North  East,  and  vertical.  The 
Residual  for  the  Vertical  is  always  zero  because  the  Vertical  is  only  used  for 
finite  differences  and  therefore  is  not  fitted.  It  can  be  seen  that  the  only 
significant  reduction  of  the  Residual  occurs  when  the  linear  effect  is  removed, 
i.  e.  ,  the  linear  fit  is  significantly  better  than  the  average.  In  order  to  verify 
these  results,  the  first  flight  was  rerun.  The  sums  of  squares  needed  in  the 
ANOVA  table  were  printed  out  with  the  altitude  and  number  of  points  preceding 
each  group.  A  sample  of  this  printout  is  shown  in  Figure  5.  Each  group  con¬ 
sists  of  10  lines.  These  correspond  to:  the  total  sums  of  squares;  that  part 
accounted  for  by  the  linear  term;  the  Dart  accounted  for  by  the  quadratic  etc., 
up  to  the  ninth.  The  criteria  for  stopping  the  fit  at  a  certain  degree  was 
removed  for  this,  allowing  each  of  the  three  coordinates  to  go  to  the  ninth 
degree.  When  the  horizontal  components  of  the  entire  flight  were  considered, 
approximately  75%  of  the  sum  of  the  squares  of  the  degree  higher  than  first 
were  small,  i.e.  ,  (did  not  explain  much  of  variation),  and  about  95%  of  the 
sums  of  the  squares  of  the  degrees  higher  than  the  second  were  not  significant. 
The  vertical  component  was  almost  never  significant  above  the  first  degree. 
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Figure  4.  Sample  Printout  of  Ninth  Degree  Method  Figure  5.  Sample  Printout  of  Modified  ANOVA  Table  for  Ninth  Degree  Method 


With  this  in  mind,  the  plots  of  the  raw  data  were  studied  again.  Almost 
without  fail,  whenever  the  degree  of  the  polynomial  above  the  second  was 
required  (by  the  tracking  error  criteria),  the  coordinate  value  plot  showed 
either  a  rough  section  of  data  or  a  relatively  smooth  section  with  one  or  tv.o 
slightly  skew  points. 

These  findings  were  the  basis  for  discarding  this  long  and  rather 
elaborate  method  in  favor  of  a  simpler  one.  The  ANOVA  table  indicated 
that  the  degree  of  the  polynomial  fit  should  probably  not  be  higher  than  two. 

1.2  SECOND-DEGREE  POLYNOMIAL  METHOD 

The  study  of  the  ninth-degree  polynomial  fit  indicated  that  in  25%  of  the 
cases  the  second  degree  polynomial  term  was  statistically  significant  in  the 
approximation.  Another  program  was  then  written  to  compute  the  wind  and 
shears  based  on  the  second  degree  fit.  Similar  expressions  were  used  for  the 
winds,  but  the  original  form  of  the  shear  was  used  because  the  component 
accelerations  were  significant  in  the  second-degree  fit.  Furthermore,  since 
the  degree  of  the  polynomial  was  reduced,  the  sampling  frequency  was  also 
reduced  to  0.2  second. 


1.2.  1  ORTHOGONAL  POLYNOMIALS  AND  DERIVATIVES 


The  orthogonal  polynomial  for  the  second-degree  fit  will  be  repeated 
here  for  ease  of  reference,  and  the  second  derivative  will  be  added.  The 
number  of  points  in  the  layer  is  n+1,  and  the  index  (i)  takes  on  all  integral 
values  from  zero  to  n.  The  Legendre  polynomials  are: 

(33a) 
(33b) 
(33c) 


(34a) 

(34b) 

(34c) 


Pn  n(l) 
0,  n 


P.  (i)  =  1  -  — 
1 ,  n  n 


P  (i)  =  1  -  —  +  j1'1.), 

2,  n  n  n  (n-1) 


The  first  derivatives  are: 


P'o.n*1*'0  P,2.n«‘>  =7T-nl^n+„-Fm 


The  second  derivatives  are: 
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P"  (i)  =  0  P".  (i)  =  C  P"  (i)  =  — t— — yr 

0,  n'  '  1,  nw  2,  n '  '  n(n-l ) 


The  first  derivative  (the  velocity)  is; 


=  AlP’l,n^  +  A2P,2,n<l> 


(35a) 

(35b) 

(35c) 


(3b) 
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The  second  derivative  (the  acceleration)  is: 


-  A2P"2,n<S> 


where  in  all  cases 


E  q.  P  (i) 

b uo  1  m’n 

lm  ”  n  ? 

E  P  (i) 

ixO  m*n 


(37) 


(12) 


1.2.2  WINDS 

In  the  second-degree  fit,  as  in  the  ninth,  we  can  compute  two  values  for 
wind.  One  is  an  average  over  the  layer,  and  the  other  is  assigned  to  a 
particular  altitude  point  in  the  layer  known  as  the  average  wind  and  point  wind 
respectively.  In  working  with  the  second-degree  polynomial,  the  values  of 
the  two  types  of  winds  are  exactly  the  same.  This  is  due  to  the  property  of 
the  Legendre  polynomial  of  the  even  degrees  being  zero  when  evaluated  at 
the  mid-point.  Thus,  since 

P,2,n<7>  "  0  <38> 


then 


alp' 


(39) 


The  quadratic  term  is  lost,  and  the  answer  is  the  same  as  if  we  had 
fitted  a  linear  equation.  However,  for  the  linear  equation  the  slope  is  the 
same  at  any  point  in  the  layer  and  is  equivalent  to  the  ratio  of  difference  in 
space  coordinate  to  the  difference  in  time,  i.e., 


slope  =  ^3  =  'ft 
r  At  q 


W 

^n/2 


(40) 


The  equations  for  the  component  winds  are: 


W 


x 


W 

y 


(41a) 

(41b) 


where  the  10  (instead  of  2)  is  necessary  to  convert  to  the  correct  units  of 
ft/sec  from  ft/0.2  sec  for  0.2  sec  spaced  data. 
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1.2.3  SHEARS 

Since  the  second  degree  was  often  found  to  be  significant,  it  was  decided 
to  use  the  following  form  for  shear: 


(9a) 

(9b) 


This  equation  implies  that  within  the  layer  there  is  a  constant  acceleration  in 
the  horizontal  components  and  a  constant  velocity  in  the  vertical  component. 
The  equations  for  shears  of  the  component  directions  in  terms  of  the  polyno¬ 
mial  coefficients  are: 


where  the  factor  5  converts  to  normal  shear  units  of  sec  .  Substitution 
for  the  Legendre  polynomials  gives; 


S 

y 


A. 

lz 


(43a) 

(43b) 


where  again  -30  is  used  in  place  of  -6  to  correct  for  the  0.2  sec  interval. 

1.2.4  ERROR  ANALYSIS  OF  WINDS  AND  SHEARS 

The  second-degree  method  does  not  compare  the  residual  to  the  tracking 
error,  as  does  the  ninth  degree.  The  error  in  the  value  of  a  parameter  was 
included  with  the  parameter  itself  in  the  output  of  :he  program  to  provide 
confidence  limits  on  the  parameters. 

1.2.4.  1  Standard  Error  of  Estimate  of  the  Coefficients 


The  standard  error  of  estimate  (cr^  )  of  the  m  th  polynomial  coefficient 
(An  )  is  given  by  m 


n 

E 

i=0 


P 

m, 


(44) 


where  cr  ^  is  the  residual  variance  of  the  m  th  degree  fit.  Therefore,  the 
squares  of  the  error  in  Aj  and  A^  are  given""By 
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(T 

e 


1 


2 


n 

£ 

i  =  0 


(7 


2 


e2 


n 

£ 

i=0 


(45a) 

(45b) 


by 


The  expression  for  the  residual  variance  of  the  first  degree  fit  is  given 


2  4  *  <qt  - 

j  L  i=o  1 


n 


q)Z  -  A/  E  P. 

1=0  1 


>]/[■' 


(46) 


and  for  the  second  degree  by 


7  f  n  7  2  n 

%  =  I  E  (qi  ‘  q)  ‘  A1  £ 

l2  [i  =  0  1  1  1  =  0 


n 


1,  n 


2  2  “  2 
(i)-A/  £  P,  J(i) 


i  =  0 


2,  n 


n -2 


(47) 


1.2.4. 2  Errors  in  Winds  and  Shears 


If  the  errors  in  the  coefficients  of  the  fitting  polynomial  are  known,  the 
errors  in  the  wind  and  shears  can  be  computed.  The  expression  for  the  wind 
is  given  by 


W 

q 


AlPVn<7> 


(48) 


and  the  error  in  the  wind 


(aw)  is  Riven  by 


a 


W 


( - )  cr 

1  n  1  A 


since  P' .  (  ^  )  =  -  2/n 

1 ,  n  c. 


(49) 


The  wind  can  then  be  reported  as  a  value  plus  or  minus  the  error  by 
using  the  following  expression: 


W 


10 

n 


,  a  1 


(50) 
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The  expression  for  shear  is  given  by 


Sq  -  q  n 


(51) 


Assuming  that  the  interaction  between  the  component  acceleration  and 
the  vertical  velocity  is  zero,  the  error  in  shear  (<r  )  is  given  by 

9 


1  2  ,  q 

=  77  ^  + 

z  n  z 


(52) 


The  square  of  the  error  in  the  component  acceleration  (a*  )  is  given  by 


\  *  [k^tt] 


2  »  2 
o-A 

2q 


(53) 


since 


’■[jA] 


2q 


The  square  of  the  error  in  the  vertical  velocity  (cr.^)  is  given  by 


2  4  —  *  2 
O'.  =  — *-  oa  . 

z  2  lz 


(54) 


n 


since  z  s 


-i  A. 

n  1  z 


Making  these  substitutions  the  error  (<r  )  becomes 

s 


(55) 


Simplifying  this  expression  gives 


•  oA 


1/2 


(5b) 
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The  shear  can  then  be  reported  as  a  value  plus  or  minus  the  error, 
using  the  following  expression: 


1.2.4.  3  A  NOVA  Table  for  Second -Degree  Fit 

The  analysis  of  variance  table  was  constructed  to  determine  the  signif¬ 
icance  of  each  term  in  the  second-degree  fit.  A  sample  of  an  ANOVA  Table 
was  given  in  Figure  2.  For  the  second-degree,  the  table  contains  only  four 
lines:  the  total  sums  of  e^uares,  that  part  due  to  the  linear  term,  that  part 
due  to  the  quadratic  term,  and  the  residual.  For  ease  of  computation,  the 
linear  and  quadratic  effect  terms  can  be  simplified  by  the  following  relation¬ 
ship. 


The  sum  of  squares  for  the  m  th  degree  is  given  by 
[  £  q  P  (i)  1  /  f"  £  P  (i)2] 

[_i = 0  1  m,n  J  [  i=0  m,n  J 


since  A 


m 


=  £  q  P  (i)  1  2  /  [  £  P  (i)2  1 

Lito  1  m*n  J  Luo  ™’n  J 


The  sum  of  squares  is  just 


A  2  [  £  P  (i)2] 

m  [  u  0  m’ n  J 


(58) 


The  A  's  are  computed  in  the  fits.  The  second  factor  is  a  constant  for 
a  given  numt&r  of  points  and  will  therefore  be  designated  as  SSI  for  the  first 
degree  and  SS2  for  the  second.  The  sums  of  squares  column  then  becomes 


(1) 

n  —  2 

E  (q:  -  q) 

i=0 

Total 

(2) 

Aj2  *  SSI 

Linea  r 

(3) 

Aj2  *  SS2 

Quadratic 

H) 

Residual  =  (1)  - 

(2)  -  (3) 
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1.2.  5  ELIMINATION  OF  STRAY  DATA  POINTS 


It  is  always  a  problem  with  a  completely  automatic  system  to  eliminate 
the  obviously  bad  data  before  it  can  contaminate  the  results  of  the  analysis. 
This  obviously  bad  data  can  be  the  result  of  a  variety  of  events  such  as 
momentary  radar  failure,  read-out  failure,  a  bad  spot  on  the  radar  tape, 
computer  error,  etc.  The  condition  usually  exists  for  only  one  data  point, 
and  the  result  is  a  value  of  a  coordinate  that  bears  little  relationship  to  the 
surrounding  points.  If  the  error  is  in  the  radar  system,  the  stray  point 
contaminates  all  of  the  coordinates;  if  it  occurs  in  computer  processing,  it 
usually  contaminates  only  one  of  the  coordinates. 

Scoggins  uses  the  common  technique  of  fitting  a  simple  function  to  the 
data  and  rejecting  a  point  if  it  falls  outside  a  tolerance,  replacing  it  by  a 
predicted  value.  As  he  points  out,  this  method  is  adequate  and  requires 
little  computer  time. 

In  our  previous  work  with  radar  tracking  data,  we  have  developed  a 
simpler  and  faster  method  of  editing  the  radar  tapes^.  It  is  based  on  the 
idea  that,  in  the  great  majority  of  cases,  there  are  only  one  or  two  con¬ 
secutively  bad  points,  and  the  surrounding  data  is  normal.  If  the  maximum 
velocity  in  all  components  is  known,  then  the  maximum  average  change  in  any 
component  from  one  point  to  the  next  can  be  computed.  This  point-to-point 
difference  is  multiplied  by  10  as  a  safety  factor  against  discarding  valid  data 
to  establish  a  maximum,  reasonable  point-to-point  change.  This  maximum 
value  is  called  the  tolerance.  In  the  editing  process,  finite  differences  are 
computed  and  compared  to  the  component  tolerances.  If  the  difference  is 
less  than  the  tolerance,  the  point  is  accepted  as  valid.  If  the  difference  is 
greater,  the  point  is  discarded,  and  the  next  point  is  read  in  and  checked  to 
see  if  the  new  differences  are  less  than  twice  the  tolerances  (the  time  interval 
has  now  doubled).  If  the  difference  is  still  larger  this  point  is  discarded;  the 
next  point  is  then  read  in,  checked  against  three  times  the  tolerance,  etc. 
until  a  "good"  point  is  found.  The  discarded  values  are  then  replaced  with 
values  computed  by  linear  interpolation  between  two  good  points.  Even  if 
only  one  coordinate  is  in  error,  the  complete  data  point  is  replaced.  If  a 
maximum  number  (50)  of  interpolated  points  is  exceeded,  the  technique 
considers  that  a  loss  of  track  had  occurred;  the  flight  is  therefore  processed 
in  separate  batches  each  having  continuity  of  track.  This  method  also  allows 
the  program  to  discard  data  in  a  time  inversion  or  to  fill  in  a  small  portion 
of  missing  data. 

While  this  method  has  less  statistical  basis  than  the  method  commonly 
employed,  it  has  proved  very  satisfactory  in  our  use  and  takes  very  little 
time  to  perform. 

1.2.6  RESULTS  AND  CONCLUSIONS  FOR  SECOND-DEGREE  METHOD 

The  basic  program  was  written  as  outlined  in  this  section.  It  contains 
the  editing  procedure  and  the  computations  for  the  wind  and  the  shears.  Two 
minor  modifications  were  made  to  the  basic  program  in  actual  operation;  one 
was  the  addition  of  the  errors  in  the  coefficients  and  the  variables  necessary 
for  the  ANOVA  Tables;  the  other  was  the  addition  of  the  errors  in  winds  and 
shears  to  the  computer  output. 
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A  printout  of  a  test  flight  with  the  modifications  of  computing  the  errors 
in  the  coefficients  and  parameters  in  the  ANOVA  Table  is  presented  in 
Figure  6.  Each  altitude  layer  presented  has  four  lines  of  output  associated 
with  it  (except  the  first  1050  feet,  which  has  only  two  lines).  The  first  two 
lines  of  each  printout  are  bracketed  and  are  the  variables  for  use  in  the 
ANOVA  Table  (these  two  lines  are  missing  for  the  first  printout).  The  last 
two  lines  contain  the  normal  printout. 

The  variables  printed  out  for  use  in  the  ANOVA  Table  are: 

1st  line  A.  ,  A,  ,  A.  ,  A,  ,  A. 

lx  2x  ly  2y  lz 

2nd  line  1)  Total  sums  of  Squares  for  x 

2)  Total  sums  of  Squares  for  y 

3)  Total  sums  of  Squares  for  x 

4)  £  P.  (i)2  =  SSI 
i=0  l'n 

5)  £  P,  (i)z  =  SS2 

i  =  0  ^,n 

The  ANOVA  Table  can  be  constructed  using  this  information  and  the 
value  of  "N"  in  the  normal  output.  The  normal  printout  consists  of  two  lines: 
the  fist  line  having  its  title  at  the  top  of  the  figure,  and  the  second  line  having 
its  title  at  the  bottom  of  the  figure. 

The  ANOVA  Tables  were  constructed  for  the  first  46  altitude  layers  for 
each  coordinate,  with  the  following  results: 

X  Y 

33  38  first  degree  only  significant 

13  6  both  first  and  second  degree  significant 

0  2  second  degree  only  significant 

This  indicates  that  since  the  second -deg ree  term  is  not  significant  in 
at  least  70%  of  the  layers,  the  second  derivative  is  not  representative  of  the 
balloon's  acceleration. 

Further  evidence  of  this  is  provided  by  examining  the  one  sigma  level  of 
the  shears.  These  va  ues  are  printed  out  as  normal  output  in  the  second 
modification  of  the  basic  program.  A  sample  of  this  output  is  presented  on 
Figure  7.  The  sigmas  of  the  shears  appear  as  the  sixth  and  seventh  data 
fields  in  the  second  card  whose  title,  as  before,  appears  at  the  bottom  of 
the  figure.  Upon  examining  the  relative  magnitude  of  the  shear  and  its  one 
sigma  level,  it  becomes  evident  that  in  some  cases  the  one  sigma  value  is 
larger  than  the  value  of  the  shear  itself.  In  fact,  in  most  levels,  the  value 
of  sigma  is  much  too  large  to  produce  confidence  in  the  values  of  shear. 


=  E  (x  -  x)Z 

=  E  (y  -  y)Z 

=  E  (z  -  t)z 
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F  gur*  6  Sample  Printout  of  Second  Degree  F.gure  7.  .-ample  Printout  of  Second  Degree  Method  with  One  Sigma  Le*e 

Method  with  Modified  ANOVA  Table  for  Wind*  and  Shear* 


The  winds  computed  by  the  second -degree  fit  agree,  for  the  most 
part,  with  the  winds  produced  by  the  ninth-degree  method,  but  the  shears 
produced  by  the  two  methods  do  not  agree  with  each  other  or  with  values 
generally  thought  to  be  realistic. 

In  general,  if  the  second  degree  is  not  significant,  the  accelerations 
predicted  by  the  fit  are  smaller  in  magnitude  than  the  noise  level  of  the  data. 
The  noise  can  be  attributed  either  to  the  actual  random  motion  of  the 
balloon  itself  or  to  the  random  tracking  errors  of  the  radar  which  will 
produce  an  apparent  random  motion  of  the  balloon.  It  is  impossible  to 
determine  which  of  these  factors  is  the  primary  cause  of  the  noise.  Thus, 
we  must  conclude  that  with  the  present  data,  only  linear  smoothing  is  valid. 

1.3  LINEAR  FIT  METHOD 

The  results  of  the  work  with  the  second-degree  polynomial  indicated 
that,  in  most  cases,  the  linear  fit  was  sufficient  to  describe  the  motion  of 
the  balloon.  Since  an  acceleration  term  could  no  longer  be  computed,  it  was 
necessary  to  go  back  to  the  finite  difference  between  winds  to  compute  the 
shear.  In  computing  the  linear  fit,  one  of  the  steps  produces  the  mean  over 
the  interval.  A  finite  difference  between  this  mean  and  the  one  from  the 
previous  layer  gives  a  type  of  average  wind  over  a  200-foot  layer.  A  shear 
can  then  also  be  computed  from  these  winds.  It  was  also  desired  to  compute 
the  one  sigma  level  of  the  linear  winds  and  include  it  in  the  printout.  The 
sampling  frequency  was  0.2  sec  for  this  degree  also. 

1.3.1  WINDS 

The  orthogonal  polynomials  for  the  first-degree  fit  were  given  in 
Section  1.2.  1 . 

In  the  linear  fit  the  point  wind  and  the  average  wind  in  the  layer,  as 
defined  in  the  section  on  the  ninth-degree  polynomial  are  identical.  In  this 
section,  we  will  introduce  a  slightly  different  definition  of  the  average  wind 
which  will  include  data  from  a  deeper  layer  and  will  be  different  from  the 
point  wind. 

1.3.  1.1  Point  Winds  and  Errors 

The  point  wind  is  the  first  derivative  evaluated  at  the  mid-point  of  the 
layer.  The  equations  for  the  component  winds,  derived  in  Section  1.2.2,  are 


W 


x 


The  errors  derived  in  Section  1.2,4  are 


(41a) 

(41b) 


(59) 
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The  point  wind  can  then  be  reported  as  a  value  plus  or  minus  the 
error  using  the  following  expression 
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1_0 
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Ai  * 


H) 
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(50) 


1.3.  1.2  Average  Winds 

When  a  linear  fit  is  performed,  the  simple  mean  of  the  data  points  in 
the  layer  is  obtained.  It  was  decided  to  compute  a  wind  from  the  finite 
differences  of  these  values  from  consecutive  layers.  Hence,  for  this  sec¬ 
tion  the  definition  of  the  average  wind  will  be: 


1+  100 


-  V/A* 


(60) 


where  *■  is  the  altitude  of  the  center  of  a  certain  layer,  and  f  +  100  is  the 
altitude  of  the  center  of  the  next  100-foot  layer.  At  is  of  course  the  time 
difference  between  q^+ioO  anc*  The  value  reported  for  the  wind  is  then 

actually  using  data  over  a  200-foot  layer.  As  one  would  expect,  this  deeper 
layer  has  a  smoothing  effect  on  the  winds  which  also  reduces  the  magnitude 
of  the  shear. 

1.3.2  SHEARS 

Since  there  are  two  distinct  types  of  winds  computed  by  this  method,  it 
was  desired  to  compute  the  shear  associated  with  each  type  wind. 

1. 3.  2.  1  Shears  from  Point  Winds 

Since  the  second  derivative  is  not  obtained  in  this  method,  the  alternate 
method  of  computing  shears  is  used: 
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where  again  the  l  is  the  altitude  of  the  mid-point  of  the  layer. 
The  component  shears  are  given  by  the  equations: 


W  -  W  W  -  W 
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(62b) 


The  point  winds  are  computed  for  the  mid-point  of  the  layer--for  example,  at 
altitudes  of  1050,  1150,  1250  feet,  etc.  The  shears  are  then  assigned  an 
altitude  exactly  between  the  two  winds  used  in  their  computation.  For  the 
example  above,  the  shears  will  be  reported  at  1100,  1200  feet,  etc. 
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1.3 .2.2  Shears  from  Average  Winds 

The  equations  for  the  shears  from  average  winds  are  identical  to  the 
shears  from  point  winds.  However,  the  altitudes  to  which  the  winds  and 
shears  are  assigned  are  reversed.  Since  the  winds  are  corr  outed  from  the 
mean  space-coordinates  assigned  to  the  mid-point  of  the  lay  r  (e.  g.  ,  1050, 
1150,  1250  etc.  ),  the  average  winds  are  assigned  to  altitudes  of  1  100,  1200 
etc.  The  shears  are  then  assigned  to  the  mid-point  of  the  wind  data  so  that 
the  altitudes  of  the  shears  are  1150,  1250  etc. 

1.  3.  3  RESULTS  AND  CONCLUSIONS  FOR  THE  FIRST  DEGREE 


The  computer  program  was  written  as  described  in  the  previous  sec¬ 
tions;  it  included  the  Technique  for  Eliminating  Stray  Data  Points  given  in 
Section  1.2.5.  Also,  a  test  run  of  the  check-out  flight  was  made.  A  sample 
of  the  result  is  presented  in  Figure  8.  The  title  at  the  top  of  the  page  is  for 
the  data  appearing  on  the  odd  lines  (those  with  50  foot  altitudes),  while  the 
title  at  the  bottom  is  for  the  data  on  the  even  lines  (those  with  100  foot 
altitudes).  The  parameters  for  the  50  foot  altitudes  are: 


Pa  ramete  r 
Alt 

Wind  East 
Wind  North 
Wind  Total 
Wind  Angle 
Average  East 
Average  North 
Sigma  East 
Sigma  North 
N 


Description 

Altitude 

Point  wind  in  East  direction 
Point  wind  in  North  direction 
Point  wind  Vector 
Angle  of  wind  Vector 

Mean  of  East  components  in  this  layer 

Mean  of  North  components  in  this  layer 

One  sigma  level  of  wind  East 

One  sigma  level  of  wind  North 

Number  of  0.2  sec  data  points  in  this  layer 


The  parameters  for  the  100  foot  altitudes  are: 


Pa  ramete  r 

- XTt - 

East  Shea r 
North  Shear 
Total  Shear 

Avg  WE 
Avg  WN 
Shear  Avg  WE 
Shear  Avg  WN 
N 


Desc  ription 
Altitude 

Shear  from  East  Point  wind 
Shear  from  North  Point  wind 
Magnitude  of  Vector  sum  of  East  and 
North  Shears 

Average  wind  in  East  Direction 
Average  wind  in  North  Direction 
Shtar  from  Avg  WrE 
Shear  from  Avg  WN 

Number  of  0.2  sec  data  points  in  this  layer 


The  point  winds  of  this  method  are  identical  to  the  point  winds  of  the 
second-degree  method,  as  they  should  be.  The  shears  from  the  point  winds 
seem  to  be  less  erratic  than  the  shears  obtained  from  the  second-degree 
method,  but  they  are  still  sufficiently  unstable  to  cause  some  questions  about 
their  reliability.  The  average  winds  are  smoother  than  the  linear  winds-* 
mainly  due  to  the  deeper  layer  over  which  the  wind  is  measured.  The  shears 
from  the  avc  rage  winds  arealso  smoother(as  expected)  but  still  somewhat  c  rratic . 
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Figvir.  ».  Semple  Printout  of  Firet  Degree  Method  Figure  9.  Semple  Printout  o(  Finite  Difference  Method 


Since  the  data  was  shown  to  be  linear  for  the  most  part,  it  was  decided 
to  compute  the  winds  from  first  finite  differences  of  the  0.  2  sec  x,  y,  z 
coordinate  data  and  to  compare  them  with  this  linear  method. 

1.4  FINITE  DIFFERENCES 

The  conclusions  reached  in  the  previous  sections  indicated  that  the  data 
were  linear  for  the  most  part  with  respect  to  time.  If  the  data  were  truly 
linear,  the  slope  with  respect  to  time  would  be  identical  to  the  first  finite 
difference.  It  is  faster  (and  therefore  cheaper)  to  compute  finite  differences 
than  to  fit  a  linear  function  and  evaluate  the  first  derivative.  The  finite 
difference  winds  and  shears  should  be  similar  to  the  linear  method. 

1.4.1  WINDS 


The  winds  in  terms  of  finite  differences  are  given  by: 
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(63b) 


since  the  data  are  equally  spaced  (0.2  sec  between  points)  equations  can  be 
rew  ritten  as: 
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where  the  "  5  "  is  to  convert  to  units  of  ft/sec. 

1.4.2  SHEARS 

The  shears  are  computed  using  the  alternate  method: 

W  -  W 
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(64b) 


(61) 


where  again  f  is  the  altitude  of  the  mid  point  of  the  layer. 

The  component  shears  are  given  by: 
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for  a  100-foot  layer. 
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1.4.3  ELIMINATION  OF  STRAY  DATA  POINTS 


The  routine  used  to  eliminate  stray  data  points  in  the  finite  difference 
method  is  identical  to  that  used  in  the  quadratic  and  linear  methods.  The 
inclusion  of  a  stray  data  point  in  a  least-squares  fit  will  cause  a  very  bad 
value.  In  the  finite  dif‘erence  me.thod,  however,  a  stray  point  will  have  no 
effect  on  the  value  of  the  winds  and  shears  provided  it  is  not  the  first  or 
last  point.  There  are  on  the  average  about  twenty-five  data  oomts  in  a 
100-foot  layer;  hence  the  probability  that  one  data  point  is  stray  is  1/25 
or  .04.  Since  there  are  normally  few  stray  points  in  the  data,  the  elimina¬ 
tion  of  stray  data  could  have  been  skipped  without  great  harm  to  the  results. 

1.4.4  RESULTS  AND  CONCLUSIONS  FOR  FINITE  DIFFERENCES 

The  computer  program  was  constructed  as  indicated  in  this  section. 

A  sample  of  the  output  is  presented  as  Figure  9.  As  before,  the  title  at 
the  top  of  the  figure  is  for  the  odd  H00  foot)  lines,  while  the  title  at  the 
bottom  is  for  the  even  (50  foot)  lines.  The  results  agree  well  with  the  linear 
method  which  indicates  that  this  cheaper  method  can  be  used  instead  of  the 
linear  method  with  little  change  in  the  results.  We  base  this  conclusion, 
however,  on  examination  of  only  a  few  soundings  by  unmodified  ROSEs. 

Such  opinion  is  therefore  subject  to  possible  future  revision. 

1.5  GENERAL  COMPARISON  OF  ALL  FITS 

In  order  to  summarize  the  results  of  the  various  type  fits,  two  portions 
of  the  test  flight  processed  by  each  method  have  been  listed  side  by  side  for 
easy  comparison.  The  two  portions  chosen  were:  l)  the  lowest  part  of  the 
flight  which  was  the  most  erratic;  2)  a  portion  near  50,000-feet  altitude  which 
was  the  smoothest.  The  east  component  was  selected  because  its  magnitude 
was  larger,  and  it  was  generally  smoother.  Tables  1  and  2  are  wind  com¬ 
parisons,  while  Tables  3  and  4  are  shear  comparisons.  Tables  l  and  3 
are  the  low  altitude,  and  Tables  2  and  4  are  the  high  altitude.  I  i  Tables  l 
and  2,  the  column  "Deg"  indicates  the  highest  degree  used  by  the  m  th  degree 
polynomial;  the  column  "ANOVA"  indicates  the  highest  significant  degree. 

The  one  sigma  level  is  included  for  the  linear  and  quadratic  winds  and  the 
quadratic  shears. 

The  values  for  the  winds  obtained  by  all  the  fits  agree  well  with  each 
other  and  with  rawinsonde  data.  Since  the  method  of  finite  difference  is  the 
quickest,  it  was  decided  to  use  this  method  in  future  work.  The  values  for 
the  shears  seem  to  depend  completely  upon  the  fitting  method  used.  No 
method  seemed  to  produce  values  which  were  consistent  from  point  to  point 
even  within  that  method. 
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TABLE  1.  Comparison  of  Winds  at  Low  Altitudes  TABLE  L.  Comparison  of  Winds  at  High  Altitudes 
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2. 


SELF- INDUCED  BALLOON  MOTION 


Wind  and  shear  values  produced  by  the  various  smoothing  techniques 
indicated  that  the  winds  were  very  erratic  at  the  bottom  of  the  flight.  This 
condition  persisted  until  about  40,  000  feet,  where  it  suddenly  vanished  and 
the  wind  became  smooth.  Concurrently  with  our  work,  Murrow  and  He^jiry 
repeated  an  experiment  first  performed  in  1958  on  RaJiosonde  balloons’. 

The  new  experiment  consisted  in  releasing  ROSE  and  other  smooth  ballr  ons 
in  still  air  and  tracking  them  with  cameras.  The  results  indicated  that  the 
spherical  balloons  experienced  self-induced  motions  similar  to  those 
experienced  by  the  ea rlie r  Radiosonde  balloons.  The  direction  of  these 
motions  is  apparently  random,  but  their  magnitude  seems  to  be  a  function 
of  the  vertical  terminal  velocity.  A  plan  view  of  six  balloon  ascents  pre¬ 
sented  in  Murrow  and  Henry's  paper  shows  that  the  balloon  behaves  very 
erratically  and  generally  unpredictably  .  When  these  results  were  made 
known,  we  plotted  several  actual  ROSE  launches  from  Eglin  AFB.  From 
concurrent  Radiosonde  releases  we  expected  a  predominantly  westerly 
wind.  The  East  vs  North  plot  showed  that  superimposed  on  the  eastward 
motion  were  various  loops  and  spirals  similar  to  those  found  by  Murrow 
and  Henry. 

The  reason  for  the  self-induced  motion  has  not  been  completely 
explained,  but  it  is  generally  thought  to  be  due  to  a  type  of  turbulence  in 
the  wake  of  the  balloon.  Some  work  has  been  done  on  this  problem  and  is 
presented  as  Appendix  2  and  Appendix  5,  but  a  satisfactory  explanation  has 
not  been  achieved. 

The  mathematical  techniques  selected  for  computing  wind  and  shear 
values  were  based  on  the  initial  assumption  that  the  balloon  motion  was 
very  nearly  the  same  as  the  motion  of  the  wind.  The  self-induced  motions 
caused  this  initial  assumption  to  be  invalid.  It  is  felt  that  the  mathematical 
techniques  do  describe  the  motion  of  the  balloon;  unfortunately,  the  values 
are  not  also  the  wind  and  shear  values. 

Two  corrective  approaches  were  taken  to  make  the  ROSE  balloons 
useful  wind  sensors.  The  first  corrective  approach  was  an  attempt  to 
measure  and  analyze  the  self-induced  motions  by  means  of  a  Power  Spectral 
Density  Technique.  The  aim  of  this  approach  was  to  develop  a  model  of  the 
self-induced  behavior  for  the  purpose  of  removing  its  effect  from  the  data. 
The  result  of  this  approach  is  presented  in  Section  3.  The  second  (and 
eventually  more  successful)  approach  was  a  physical  modification  of  the 
balloon  itself  to  control  the  wake  turbulence. 
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3. 


POWER  SPECTRAL  ANALYSIS  OF  BALLOON  MOTION 


The  East  vs  North  plots  of  the  ROSE  data  indicated  that  the  balloon 
motion  was  generally  unpredictable  in  a  small  interval.  However,  when  a 
plot  of  an  entire  flight  was  viewed  there  appeared  to  be  a  type  of  fairly 
regular  oscillation.  The  technique  commonly  used  to  determine  oscillatory 
changes  in  data  is  the  Power  Spectral  Density  Analysis.  The  basis  for  this 
technique  is  given  in  the  standard  work  on  time  series  analysis  by  Blackman 
and  Tukey  .  An  excellent  Time  Series  Analysis  program  has  been  prepared 
by  Healy  and  Bogerr’  at  Bell  Telephone  Laboratories  and  is  available  through 
SHARE'.  The  final  output  of  the  program  is  a  table  of  Power  at  specific 
frequencies  vs  the  frequencies.  A  spike  in  the  curve  indicates  a  predomi¬ 
nance  of  that  frequency  in  the  data. 

3.  1  TIME  SERIES  ANALYSIS  COMPUTER  PROGRAM 

Because  of  the  sizable  programming  effort  required  to  prepare  a  Time 
Series  Analysis  for  a  computer,  a  search  of  the  available  program  libraries 
was  made.  From  the  several  good  methods  available,  ’  °  the  one  by 

Healy  and  Bogertb  was  selected  because  their  approach  seemed  easiest  to 
tailor  to  our  particular  problem.  The  time  series  analysis  involves  a  large 
number  of  computations  on  a  fairly  large  quantity  of  data.  Healy  and  Bogert 
broke  down  the  computations  into  separate  numerical  processes  and  wrote  one 
FORTRAN  subroutine  for  each  process.  The  method  was  then  to  use  each 
subroutine  (with  the  necessary  linkages)  as  a  high  le/el  command.  For 
example,  to  read  the  data  into  the  machine  the  main  program  uses  merely 
CALL  READIN  (A,  NA,  ANAME)  to  transfer  the  data  to  the  array  "A"  where 
it  is  ready  for  computation.  The  main  calling  sequence  gives  the  overall  con¬ 
cept  of  the  computation  steps  without  regard  to  all  of  the  FORTRAN  mechanics 
necessary  for  the  actual  computer  operation.  This  makes  modification  of  the 
program  very  easy.  One  merely  w  rites  and  debugs  a  new  independent  sub¬ 
routine  including  only  the  "CALL”  line  in  the  main  program. 

The  original  version  of  the  program  contained  17  subroutines  which  may 
be  found  in  Reference  6.  For  our  analysis  we  decided  to  use  just  5  of  17 
subroutines.  These  5  subroutines  are  presented  as  Appendix  3  with  descrip¬ 
tions  of  their  individual  processes. 

3.  1.  1  MODIFICATIONS  TO  THE  HEALY-BOGERT  PROGRAM 

There  were  two  major  changes  made  to  the  basic  program:  the  con¬ 
version  of  the  compiler  language  and  an  increase  in  size.  The  original 
compiler  language  used  was  FORTRAN  II,  and  all  of  the  input,  output,  and 
computations  were  coded  in  FORTRAN.  The  program  was  converted  to 
FORTRAN  IV  and  operated  under  the  IBSYS  system  in  order  to  use  the 
techniques  provided  by  this  more  advanced  system.  The  original  program 
was  limited  to  a  maximum  of  1500  points  per  analysis.  By  using  fewer  sub¬ 
routines,  we  were  able  to  increase  this  maximum  to  5000  points  per  analysis. 
The  other  minor  changes  were  the  addition  of  different  t rend- removing  tech¬ 
niques  and  the  addition  of  pre- whitening  and  post-darkening  subroutines.  The 
input  subroutine  READIN  was  also  changed  to  conform  to  our  data.  The 
additional  subroutines  are  described  in  Appendix  4. 
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3.1.2  PROCEDURE  TO  TEST  PROGRAM  USING  SIMULATED  DATA 


The  programs  received  through  SHARE  are  normally  well  documented 
and  perform  satisfactorily;  therefore  the  program  was  tested  by  simulating 
data  with  known  frequencies  and  processing  it  by  using  the  PSD  program. 

The  manufactured  data  were  a  linear  function  with  frequencies  superimposed 
on  it.  The  frequencies  were-  1  /187,  1  /131,  1  /43,  1  /41,  1  /37,  1  /17,  1  /7, 
and  1/(1  +V3)  cycles  per  second.  The  main  calling  sequence  used  the 
following  subroutines:  READIN,  DETRND  (over  the  whole  data  set), 

AUTCOV  (with  100  lags),  FOURTR  (using  the  hanned  cosine  transforms), 
and  OUTPUT.  The  program  performed  very  well  in  reporting  all  of  the 
frequencies  we  had  manufactured.  The  output  from  the  test  was  plotted  on 
rectangular  scales  (presented  as  Figure  10)  and  on  semi-log  scales 
(presented  as  Figure  11). 

The  ROSE  data  seemed  to  be  essentially  quadratic  rather  than  linear; 
hence  the  same  eight  frequencies  were  superimposed  on  a  parabola,  and  the 
PSD  program  was  re-run  with  the  QTRND  subroutine  (Appendix  4).  All  of 
the  frequencies  were  reported  again,  and  the  shapes  of  the  P  (<*>)  vs  $  curves 
were  identical  to  those  shown  in  Figures  10  and  11. 

Since  the  program  seemed  to  perform  well  on  "clean"  data,  it  was 
next  decided  to  add  to  the  quadratic  data  random  noise  of  approximately  10% 
of  the  amplitude  of  the  data.  The  pre -whitening  (WHITN)  and  post-darkening 
(PDRKN)  subroutines  (Appendix  4)  were  then  added  in  an  attempt  to  remove 
the  noise.  The  output  again  was  almost  identical  to  that  shown  in  Figures  10 
and  11. 

The  next  step  was  to  see  if  a  method  simpler  than  the  least  squares 
linear,  or  quadratic  methods  could  remove  the  general  trend  and  still  pro¬ 
duce  results  which  would  be  satisfactory.  The  first  method  (TRNDR)  was  to 
remove  a  line  determined  by  the  first  and  last  points  of  a  five-second  section. 
Then  by  using  the  last  point  of  a  section  as  the  first  point  of  the  next  section, 
continuity  could  be  produced  through  all  of  the  data.  The  results  were  again 
almost  identical  to  those  shown  in  Figures  10  and  11. 

It  was  thought  that  perhaps  a  moving  average  would  better  remove  the 
general  trend  from  the  data  than  a  series  of  straight  lines  as  in  TRNDR. 

The  subroutine  MOVAVG  was  written  and  used  with  both  51-  and  81-point 
averages.  The  number  of  lags  was  then  reduced  to  50--first  because  the 
number  of  points  was  reduced  and  secondly,  in  order  to  determine  how  the 
lower  resolution  would  aficct  the  shape  and  placement  of  the  spikes.  The 
result  of  the  51 -point  average  is  shown  as  Figure  12.  This  result  is  identical 
to  the  output  from  the  81 -point  average  method.  It  will  be  noticed  that  the 
peaks  have  broadened  slightly  but  are  still  very  satisfactory  in  returning  the 
frequencies  known  to  be  present  in  the  data. 

The  results  of  these  tests  indicated  that  the  program  was  performing 
satisfactorily  and  could  be  applied  to  actual  data. 
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Figure  13.  Azimuth  ve  Time  (or  Flight  16021  23,000  to  33,000  Feet  Altitude  Leyer 
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3.2  ANALYSIS  OF  ROSE  DATA  BY  PSD  PROGRAM 


The  vertical  motion  of  a  ROSE  balloon  is  very  smooth  and  does  not 
exhibit  any  of  the  oscillatory  motion  found  in  the  horizontal  plane.  Because 
of  this  the  Azimuth  angle  of  the  radar  was  selected  for  analysis  rather  than 
the  Elevation  angle  or  Range.  The  program  was  set  up  for  a  computer  run 
so  that  a  ROSE  packed  binary  tape  was  the  input,  and  the  PSD  table  with 
the  history  and  identification  of  the  run  was  the  output. 

3.2.1  STANDARD  ROSE  BALLOONS 

The  standard  ROSE  balloon  results  were  achieved  in  two  steps.  The 
first  was  a  control  study  using  only  one  flight  (16021)  chosen  at  random.  The 
second  step  vas  to  apply  the  techniques  developed  in  the  first  step  to  a 
scries  of  standard  balloons  flown  at  nearly  the  same  time  as  some  of  the 
modified  balloons. 

3.  2.  1 . 1  Control  Study  of  Flight  1602  1 

In  order  to  obtain  a  better  understanding  of  the  actual  ROSE  data,  a 
plot  of  Azimuth  vs  Time  for  flight  16021  was  made.  A  portion,  from 
23  thousand  to  33  thousand  feet  altitude,  is  presented  as  Figure  13.  This 
portion  of  400  seconds  is  essentially  a  quadratic  with  the  familiar  oscilla¬ 
tions  superimposed  on  it. 


The  data  were  processed  four  times,  each  using  a  different  trend- 
removing  technique.  The  main  program  used  the  following  calling  sequence: 


1.  READIN  (400  points) 

2.  One  of  the  trend-removing  routines 

3.  WHITN  (ALFA1  =  -  1/2,  ALFA2  =  +1) 

4.  AUTCOV  (50  lags) 

5.  FOURTR  (Hanned  cosine  transform) 

6.  PDRKN  (ALFAl  =  -  1/2,  ALFA2  =  +  1) 

The  four  trend- removing  routines  used  were: 

1.  QTRND  over  the  whole  interval 

2.  TRNDR  over  50  point  lines 

3.  MOVAVG  over  51  points 

4.  MOVAVG  over  81  points 

The  results  of  all  four  runs  were  almost  identical,  and  a  representative 
plot  of  P(<$>)  vs  (*>  is  shown  in  Figure  14.  The  plot  indicates  that  in  this  por¬ 
tion  of  the  flight  there  was  a  much  greater  predominance  of  low  frequencies 
than  high  frequencies,  and  as  frequency  increases  the  power  diminishes 
approximately  as  the  log  of  the  frequency.  It  is  interesting  to  note  that  the 
very  slight  peak  at  a  frequency  of  0.  1  cps  repeats  itself  in  all  four  attempts. 
However,  it  is  believed  that  the  peak  does  not  rise  high  enough  above  its 
surroundings  to  enable  one  to  say  that  there  is  a  predominant  oscillation  at 
a  frequency  of  0.  1  cps. 
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Io  determine  the  effect  of  the  oscillating  azimuth  angle  on  the 
horizontal  components,  two  analyses  were  run  (each  with  a  different  trend 
remover)  on  each  of  the  East  and  North  components  for  the  same  altitude 
band.  Although  the  balloon  was  launched  near  the  radar  site,  a  strong 
westerly  at  the  time  of  this  flight  had  pushed  the  balloon  (at  an  altitude  of 
23,000  feet)  much  farther  to  the  east  but  at  about  the  same  latitude.  This 
resulted  in  the  North  component  being  more  noticeably  oscillatory  than  the 
East  component.  The  different  t rend- removing  routines  again  had  no 
visible  effect  in  producing  a  power  curve  of  different  shape.  A  typical 
result  for  the  North  component  is  presented  as  Figure  15;  a  result  for  the 
East  component  is  shown  as  Figure  16.  The  slight  increase  in  power  is 
again  noticed  in  the  North  component  at  0.  1  cps;  because  of  the  location  of 
the  balloon  with  respect  to  the  radar,  such  increase  is  not  perceived  in  the 
East  component. 

The  next  step  was  to  work  with  various  parts  of  the  flight  analyzing 
only  the  azimuth  angle.  First,  the  entire  flight  was  analyzed  (0  to  63 ,000  feet). 
The  result  was  a  very  noisy  linear  function  on  semi-log  scale  similar  to  those 
shown  in  Figures  14  and  17.  This  function  was  not  only  noisier  than  those 
depicted  in  Figures  14  and  17,  it  also  lacked  the  prominence  at  0.  I  cps. 
Analysis  of  the  first  half  of  the  flight  (0  to  33,  000  feet)  gave  the  same  result 
as  the  analysis  of  the  whole  flight.  When  the  Azimuth  vs  Time  plot  of  the 
entire  flight  was  studied,  it  was  noticed  that  at  about  42,000  feet  the 
oscillations  vanished  (which  coincided  with  the  smoothing  of  the  winds 
mentioned  earlier).  An  analysis  of  the  last  portion  of  the  flight  (43,  000  to 
63,  000  feet)  was  made.  The  results,  presented  as  Figure  18,  indicate  the 
presence  of  very  low  frequency  (below  0.03  cps)  oscillations  and  essentially 
white  noise  at  all  frequencies  above  .this  - -a  result  quite  different  from  any 
obtained  in  previous  analyses.  This  agrees  with  the  observed  fact  that  the 
winds  suddenly  became  smooth  at  42,000  feet.  In  order  to  check  this  result, 
the  portion  of  the  flight  from  34,000  feet  to  45,000  feet  was  analyzed.  The 
result  is  presented  as  Figure  17  in  which  the  whole  first  portion  of  the  flight 
seems  to  confirm  that  the  oscillations  really  do  vanish  above  42,000  feet  for 
the  standard  ROSE  balloon. 

3.  2.  1. 2  Series  of  Standard  Flights 

The  results  of  the  control  study  indicated  that  Azimuth  should  be 
analyzed  and  that  the  flight  should  be  broken  into  two  sections:  1)  the  portion 
containing  the  oscillations  (the  first  of  the  lower  half  of  the  flight);  2)  the  part 
that  is  smoother  (the  second  or  upper  half  of  the  flight).  The  flights  of 
several  standard  balloons  were  chosen  for  analysis  of  the  azimuth  both  in 
the  two  altitude  bands  and  over  the  entire  flight.  The  flight  numbers  of  the 
balloons  chosen  were  16053,  16074,  16083,  16114,  and  16116.  The  output 
is  presented  as  Figures  19,  20,  21,  22,  and  23  respectively.  The  analysis 
of  the  entire  flight  is  not  shown  because  it  is  very  similar  to  the  lower-half 
analysis.  It  can  be  seen  that  while  the  general  trend  of  the  upper  and  lower 
portions  are  similar,  the  upper  portions  have  none  of  the  small  peaks  found 
in  the  lower  portions.  These  small  peaks  are  not  sufficiently  prominent  to 
give  assurance  that  these  frequencies  are  always  present  in  the  data.  If  the 
frequencies  are  not  always  present,  attempting  to  remove  them  could  cause 
even  more  erratic  data. 
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Figure  18.  Results  of  PSD  Analysis  of  Azimuth  of  Flight  16021  43,000  to  63,000  Feet  Altitude 
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3.2.2  MODIFIED  BALLOONS 


The  second  approach  to  eliminating  self- induced  balloon  motion  was  to 
physically  modify  the  balloons.  These  modifications  vail  be  discussed  in 
Section  4;  for  now  it  will  be  sufficient  to  indicate  that  there  were  basically 
three  types  of  modified  balloons:  tails,  fences,  and  roughening- -plus  two 
flights  by  a  one-meter  sphere  with  corner  reflectors  known  as  the  ROBIN 
Balloon.  Power  spectral  density  analyses  were  performed  on  the  two  halves 
of  at  least  one  flight  of  each  type  of  modified  balloon.  The  P(4>)  vs  <j>  plots 
are  presented  for  these  representative  flights: 


Flight  No. 

Modification  Type 

Figure 

16055 

belts 

24 

16073 

roughening 

25 

16081 

belts 

26 

16113 

tails 

27 

16115 

tails 

28 

Upper  portions  of  flights  of  modified  balloons  and  standard  balloons  are 
very  similar,  but  the  similarity  between  the  lower  half-portions  of  flights  of 
these  balloons  varies  with  the  type  of  modification  used.  In  all  cases  at 
least  some  of  the  smaller  peaks  were  removed  and  in  general  the  modified 
balloons  resulted  in  a  smoother  PSD  result. 

. 

3.2.3  MODIFIED  ROBIN  BALLOONS 

In  addition  to  the  variously  modified  ROSE  balloons  two  smooth,  one- 
meter  ROBIN  balloons  were  also  tested.  The  first  had  a  two-inch  inflation 
valve  and  a  five-mb  pressure  -  relief  valve.  The  second  had  a  two-inch 
combination  inflation  and  relief  valve.  The  flight  was  divided  into  two 
portions  with  the  division  at  23,000  feet.  The  results  of  the  PSD  analysis 
of  the  upper  portion  is  very  similar  to  the  flight  pattern  shown  in  Figure  17. 
Except  for  the  very  low  frequency  peak,  the  data  contained  the  same  power 
at  all  frequencies.  The  lower  portion  is  similar  to  the  lower  portions  of  the 
other  flights  with  one  distinction.  There  is  a  broad  peak  in  the  power  curve 
from  about  0,22  to  0.245  cps  which  rises  high  enough  above  its  surrounding 
to  give  it  a  degree  of  credibility.  The  rise  rate  of  the  ROBIN  is  about 
15  ft/sec  in  this  region;  therefore  the  oscillation  covers  60  to  70  feet  of 
altitude.  However,  the  presence  of  the  large  amount  of  power  at  lower 
frequencies  (as  seen  in  Figure  29)  makes  this  discovery  impractical  to 
apply  in  a  data  reduction  technique. 

3.3  GENERAL  CONCLUSIONS 

The  purpose  of  this  approach  was  to  find  a  predictable  periodicity  which 
could  be  eliminated  from  the  data  ii  order  to  produce  reliable  winds  and 
shears.  We  have  not  been  able  to  find  (with  any  degree  of  certainty)  such  a 
predictable  periodicity  since  the  motion  of  the  balloon  apparently  contains 
no  prominent  frequencies. 
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4. 


PHYSICAL  MODIFICATION  TO  THE  STANDARD  ROSE  BALLOONS 


The  second  approach  to  making  the  ROSE  a  more  reliable  wind  sensor 
was  to  make  various  physical  modifications  to  the  shape  of  standard  ROSE 
balloons  in  an  attempt  to  control  the  wake  separation. 

The  standard  ROSE,  as  described  previously,  is  a  smooth  sphere. 

A  photograph  of  a  standard  ROSE  taken  immediately  after  launch  is  shown 
in  Figure  30.  The  modifications  were  of  three  basic  types:  tails,  orthogonal 
fences,  and  over-all  roughening.  The  tail  modifications  were  devices  similar 
to  those  used  with  ordinary  kites  to  improve  the  stability  without  changing  the 
shape  of  the  standard  ROSE  balloon.  The  orthogonal  fence  modifications 
(shown  in  Figure  31)  consists  of  two  orthogonal  great  circle  strips  of  poly¬ 
film  with  lightweight  drinking  cups  attached  with  their  open  end  toward  the 
balloon.  The  over-all  roughening  modification  (shown  in  Figure  32)  is 
popularly  known  as  the  "Jimsphere"  and  is  generally  attributed  to  Jim 
Scoggins.  It  consists  of  lightweight  cups  attached  to  the  balloon  in  a  fairly 
symmetrical  design  similar  to  a  golf  ball  that  had  been  turned  inside  out^. 

In  addition  to  these  modified  ROSE  balloons  there  were  also  two 
ROBIN  balloons.  The  ROBIN  balloon  is  a  one  meter  diameter  smooth  sphere, 
as  shown  in  Figure  33.  They  were  modified  from  standard  configuration  by 
the  addition  of  a  pressure-relief  valve  so  they  could  rise  higher  without 
bursting. 

Each  modified  balloon  was  flown  and  tracked  in  a  program  testing 
the  effectiveness  of  the  modification.  The  testing  program  is  described  in 
Section  5. 
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5. 


TESTING  OF  MODIFIED  ROSE  BALLOONS 


In  order  to  evaluate  the  effects  of  the  ROSE  modifications,  a  series  of 
89  test  flights  were  flown  at  Eglin  AFB  where  a  standard  and  a  modified 
ROSE  balloon  were  flown  at  about  the  same  time.  A  comparison  involving 
the  variances  of  the  winds  was  made. 

5.  1  COLLECTION  OF  DATA 

Each  flight  of  a  balloon  was  tracked  with  the  AN/FPS-16  radar,  and 
the  data  was  processed  by  Eglin’s  standard  xyz  reduction  program.  When 
a  group  of  processed  flights  was  completed,  they  were  stored  serially  on 
a  binary  magnetic  tape.  The  binary  tape  was  then  forwarded  to  the 
University  of  Dayton  Research  Institute  together  with  a  BCD  listing  of 
every  tenth  pcint  (in  order  to  reduce  the  volume  of  the  paper).  An  inven¬ 
tory  of  this  test  program  is  presented  as  Table  5  and  includes  all  the 
pertinent  information  concerning  the  flights. 

5.2  COMPUTATION  OF  WINDS  AND  THEIR  VARIANCES 

The  conclusions  reached  in  the  studies  of  the  various  fits  indicated  that 
the  finite-difference  method  of  computing  wind  is  normally  adequate.  Finite- 
difference  Winds  over  100-foot  layers  of  altitude  were  computed  for  each 
flight  of  the  series.  When  the  mean  winds  over  the  1000-foot  layer  were 
computed,  the  variances  of  both  components  and  the  wind  vector  over  the 
same  1000-foot  layer  were  also  computed.  This  was  done  in  order  to  quanti¬ 
fy  the  effect  of  the  modification.  A  sample  of  this  type  of  printout  is  pre¬ 
sented  as  Figure  34.  A  plot  of  the  Wind  Vector  vs  Altitude  was  drawn  for 
each  flight.  Samples  of  portions  of  those  plots  are  shown  as  Figures  35  to  38. 
The  large  dot  shown  in  these  figures  represents  the  average  over  1000  feet 
altitudes  while  the  x'ed,  dashed  line  is  the  value  reported  by  standard 
rawinsonde  balloon  launched  at  about  the  same  time  as  the  test  balloon. 

Figure  35  describes  the  wind  reported  by  a  standard  ROSE  balloon  at  the 
altitude  at  which  oscillations  disappear.  Figures  36  to  38  are  plots  of  the 
variously  modified  ROSE  balloons  which  present  a  general  view  of  the  effect 
of  the  modifications  on  the  computation  of  the  Wind  Vector. 

5.  3  COMPARISON  OF  VARIANCES 

A  representative  variance  of  the  wind  vector  over  a  10000-foot  altitude 
layer  was  computed  using  the  following  equation: 

2  El*'"',/  + 

"■  *  - n -  •  <«) 

where  E<rWr^  and  £  o'Wjyj^  are  the  sums  of  the  individual  1000-foot  vanances 
over  each  10000-foot  layer.  These  variances  were  computed  for  four  layers 
from  0  to  40000-foot  altitude  and  are  presented  as  Table  6.  The  comparison 
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of  variances  was  performed  by  computing  the  Ratio  (IT)  given  by  the  follow¬ 
ing  equation: 


2  2 
<r  -  a 

—  _  Std  modified 
R  =  - 2 - 


(66) 


Std 


where  (r^_  .  and  <r^  .are  the  variances  over  10000-foot  altitude  layer 

Std  modified  ' 

of  a  standard  ROSE  and  a  modified  ROSE  respectively,  flown  at  nearly  the 
same  time.  The  ratio  for  each  of  the  four  altitude  layers  for  each  modified 
balloon  is  presented  as  Table  7.  The  Ratio  (IT)  is  a  type  of  percent-reduction 
equation,  and  the  values  of  R  for  a  modified  balloon  can  be  interpreted  as 
follows:  a  small  value  for  R  indicates  that  there  is  very  little  difference 
between  a  standard  balloon  and  this  modified  balloon;  a  large  value  for  7 
indicates  that_there  was  a  significant  change  due  to  the  modification.  When 
the  value  of  R  approaches  1.0  the  variance  of  the  modified  balloon  approaches 
zero--which  is  just  about  as  unrealistic  as  the  winds  reported  by  the  unmodi¬ 
fied  ROSE  balloon.  It  was  determined  that  a  0.  7  to  0.  9  reduction  in  variance 
would  indicate  a  successful  modification.  Some  of  the  more  successful 
modifications  would  then  be  applied  to  the  remaining  standard  ROSE  balloons 
in  the  Air  Force  stocks  so  they  could  be  used  to  produce  reliable  wind  values. 
Care  must  be  used  in  the  selection  of  the  modification  because  the  balloon's 
wind  response  can  be  adversely  affected.  ^ 
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Figure  34.  Sample  Output  of  Program  to  Compute  Winds  for  the  Evaluation 
of  the  Modified  ROSE's 
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Figure  35.  Wind  Vector  ve  Altitude  Flight  16021;  40,  000  to  50.  000  Feet 
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Figure  36.  Wind  Vector  ve  Altitude  Flight  16055;  ifl.OOO  to  38.000  Feet 
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Figure  37,  Wind  Vector  v«  Altitude  Flight  16073;  22.  700  to  32.  700  Feet 
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Figure  38.  Wind  Vector  ve  Altitude  Flight  16C8 1  23,600  to  33,600  Feet 
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TABLE  7.  Comparison  of  Variances  for  Modified  ROSE  Balloons 


Modified  Std 


All.  Layer 

0-10, 000 


Alt.  leaver 
1  I -20,  000 


'It.  Laser 
2  I  -  JO,  000 


Ail  Laser 
Jl  -40,  000 


16002 

16001 

0.  500 

0.  760 

0.  569 

0.  074 

16004 

16003 

0.  763 

0.  809 

0.  526 

0;  32  8 

1601  1 

16012 

N.  G. 

0.  287 

0.  406 

0.  1  79 

I601J 

16015 

0.  1  30 

0.  148 

0.  362 

0.  095 

1  602  J 

16021 

0.  806 

0  904 

0.  926 

0.  64  4 

1  603 1 

16035 

0.  919 

0.  935 

0.  907 

N'  G 

16032 

16037 

0,  717 

0.  606 

0  710 

0.  397 

1603  J 

1  60J7 

0,  878 

0.  897 

0  892 

0.  2  1  1 

I60J4 

160J9 

N.  G. 

N.  G. 

0.  ,10 

0.  426 

1  6038 

16039 

0.  507 

0.  274 

N.  G. 

N  G 

16042 

1604  1 

N,  G. 

N.  G. 

0.  638 

0.  1  51 

1604  3 

16041 

0.  736 

0.  639 

0.  671 

0.  389 

1  6045 

16044 

N.  G. 

0  535 

0.  377 

N.  G. 

16047 

16046 

0.  329 

0.  555 

0.  686 

0.  856 

16052 

16051 

0  256 

0.  216 

0.  1  50 

0.  64  9 

1' 054 

16053 

N.  G 

0.  452 

0  501 

0.  942 

16055 

I605J 

0.  7J5 

0.  705 

0.  768 

0.  657 

16036 

16057 

0.  356 

0.  248 

N.  G. 

5J  C, 

16061 

1  6  0  b  6 

0.  894 

0.  882 

0  814 

N.  G. 

16062 

16066 

0.  559 

0.  428 

0.  401 

N.  G. 

I606J 

16067 

0.  429 

0.  501 

0.  493 

0.  421 

16064 

16068 

0  659 

0.  778 

0.  708 

0.  270 

16065 

16068 

0.  1  18 

N.  G. 

0.119 

N  G. 

16071 

16072 

0  191 

0.  201 

0.  155 

N.  G 

1  607  J 

16074 

0.  690 

0.  618 

0.  599 

0  9b  1 

16075 

16074 

0.  955 

0.  860 

0  885 

0.  988 

16081 

16083 

0.  475 

0.  450 

0.  586 

0.  274 

16082 

1608  3 

0.  874 

0.  871 

0.  86  3 

0.  889 

16085 

16084 

0.  759 

0.  773 

0.  670 

0.  502 

16091 

16092 

0.  625 

0.  670 

0.  596 

0.  5J6 

1  609  J 

16004 

N.  G. 

0.  218 

N.  G. 

0  4  69 

16096 

16091, 

0  514 

0.  745 

N.  G. 

S'  Ci. 

16101 

16102 

0.  6  30 

0.  607 

0.  596 

0.  536 

16103 

16104 

N.  G. 

0.  218 

N  G 

N  G 

16106 

16  105 

N.  G. 

N.  (,. 

N.G. 

0  9|J 

161  1  1 

16112 

N  G 

0.  91  3 

0.  922 

0  842 

16113 

1  6  1  1  4 

0.  489 

0.  549 

0.  548 

0.  554 

16115 

16166 

0.  519 

0  916 

0.  55} 

V  G. 

16)2  J 

16124 

N  G. 

0.  221 

0.  838 

N.  G 

1612  5 

16126 

N  G 

N.  G. 

N.  G. 

0.  616 

16120 

161  JO 

N  G. 

N  G. 

N.  G. 

0  934 

161  32 

16131 

N  G 

0,  80} 

N.  G. 

N'.  G. 

16134 

161  3  3 

0.  8  35 

0,  66  5 

0.  591 

0.  931 

6. 


SMALL  LIGHTW  EIGHT  ROSE  BALLOONS 


The  two  appioaches  initially  suggested  to  make  the  ROSE  system 
reliable  wind  sensors  met  with  limited  success.  A  new  approach  which 
involved  designing  an  entirely  new  balloon  sensor  was  suggested  by  Reid. 
Experiments  tend  to  confirm  that  the  self -indue  .id  balloon  motion  can  be 
related  to  the  Reynolds  Number,  a  unitless  ratio  of  inertial  to  viscous  forces. 
When  the  Reynolds  Number  is  greater  than  250000  the  balloon  experiences 
considerable  instability,  but  when  the  Reynolds  Number  is  less  than  250000 
the  instability  is  negligible.  The  value  of  250000  will  be  known  as  the 
critical  Reynolds  Number.  The  change  from  supercritical  to  subcritical 
occurs  at  about  42,  000  feet  altitude  for  the  standard  ROSE  balloon.  Reynolds 
Number  for  a  sphere  is  given  by 


Re  =  (67) 

|  v|  is  the  velocity  with  respect  to  the  surrounding  air  mass,  p  is  the  density 
of  the  air  mass,  D  the  diameter,  and  p  the  coefficient  of  viscosity.  Since 
the  density  or  coefficient  of  viscosity  of  the  air  mass  is  obviously  uncontrol¬ 
lable,  Reynolds  Number  must  be  affected  by  changing  the  velocity  or  the 
diameter  of  the  balloon.  The  velocity  could  be  decreased  by  adding  weight 
to  the  balloon,  but  this  is  not  practical  because  it  also  increases  the  lag 
distance  and  the  length  of  time  required  for  a  sounding.  Decreasing  the 
diameter  of  the  balloon  lowers  the  maximum  altitude  the  balloon  can  attain 
but  improves  its  wind  response.  Reid  has  shown  that  a  balloon  having  a 
diameter  of  approximately  1  meter  and  weighing  less  than  100  gms  should 
produce  detailed  winds  from  altitudes  of  10000  to  50000  feet^®.  A  series 
of  12  such  balloons,  designated  as  GT  balloons,  were  fabricated  for  testing 
in  comparative  flight  tests  with  standard  two-meter  ROSE  balloons. 

6.  1  COMPUTATION  OF  WINDS,  VARIANCES  AND 
REYNOLDS  NUMBERS  OF  GT  BALLOONS 

The  program  to  test  the  GI  balloons  was  identical  to  the  one  used  to 
test  the  modified  ROSE  balloons.  A  list  of  the  GT  balloons  and  the  standard 
ROSEs  flown  in  the  test  series  is  presented  as  Table  8.  Wind  values  were 
computed  by  the  finite-difference  method  for  100-foot  layers  of  atmosphere. 
Average  winds  for  the  1000-foot  layers  with  their  variances  were  then  com¬ 
puted. 


In  addition  to  the  normal  output  (Figure  34),  two  extra  parameters  Rise 
Rate  and  Reynolds  Number,  were  added  to  the  printout  shown  in  Figure  39. 
Rise  Rate  was  computed  by  finite -differences ,  and  Reynolds  Number  was 
computed  from  the  following  equation: 

(Rise  Rate)  •  p  •  D 

Re  =  - F - 12 - •  <b8> 

^  m 
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uhere  p  and  p.  are  the  density  and  the  coefficient  of  viscosity  respec¬ 
tively.  T?iey  were  obtained  from  the  1962  Model  Atmosphere*  . 

In  addition,  four  flights  (16141,  16142,  16143,  and  16144)  were  recom¬ 
puted  using  the  linear  fit  method  of  computation  described  in  Section  1.3. 
Values  obtained  by  these  various  methods  were  then  compared. 

6.2  COMPARISON  OF  VARIANCES 

The  method  used  to  compare  the  variances  for  the  GT  balloons  was 
identical  (computing  the  U  values)  to  that  used  for  the  modified  ROSE 
balloons.  The  values  for  the  ratios  are  presented  as  Table  9.  The  values 
for  the  ratios  for  flights  16141  to  16144  are  presented  as  Table  10. 

The  tables  show  a  generally  high  reduction  in  the  wind  variance.  Based 
on  these  values  and  their  own  independent  calculations,  the  scientists  at 
AFCRL  have  made  a  tentative  selection  of  a  new  standard  ROSE  balloon.  It 
was  the  balloon  designated  as  GT59  in  the  test  series  and  is  a  0.25-mil, 
40-inch  balloon  with  a  standard  pop-out  valve. 
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-11.99 

15.89 

12.52 

12.51 

205460. 

16147 

5503. 

-17P6. 

-2054. 

7.  1 
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Figure  39.  Sample  Output  of  Wind  Computation  for  Standard  and  Modified 
Test  Series 
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TABLE  8.  Flight*  Conducted  In  Comparison  Te*t*  of  GT  Balloon* 


Flight*  Conducted  in  Comparison  Test  of  GT  Daiioons 


Flight 

Balloon 

Diameter 

Inches 

Gauge 

mil 

Weight  Valve 
gms 

Time  and 

Date  of  Launch 

16141 

GT  -62 

40 

.  25 

48.  1  Weightless 

1242 

26  May  »>4 

16142 

GT -19 

40 

.  21 

73.  1  Std 

1429 

26  May  64 

1 6  i  4  3 

13 

78  1 

.  50 

310  Std 

1742 

26  May  64 

16144 

30 

78.  1 

.  50 

350  Std 

1406 

26  May  64 

16141 

29 

78.  1 

.  50 

310  Std 

1128 

26  May  64 

16146 

GT -16 

36 

.  21 

29.  3  Weightless 

1629 

26  May  64 

16147 

GT  - 12 

36 

.  31 

74  9  Std 

i  608 

26  May  64 

16111 

GT-61 

40 

.  10 

7  3.7  Weightle* s 

1237 

27  May  64 

16112 

41 

78.  1 

.  10 

310  Std 

1428 

27  May  64 

1611) 

19 

78.  1 

.  50 

310  Std 

11)2 

27  May  M 

16114 

GT  - 18 

40 

.  31 

84.  8  Std 

1619 

27  May  64 

16111 

GT  -11 

36 

.  50 

62.  3  Weightless 

1411 

27  May  64 

16:  16 

GT  - 1 3 

36 

.25 

67.  2  Std 

1926 

27  May  64 

16117 

48 

78.  1 

.  10 

310  Std 

1814 

27  May  64 

16161 

GT  -17 

40 

.  10 

100.  6  Std 

1  346 

28  May  64 

16162 

2) 

78.  1 

.  10 

3  50  Std 

1211 

28  May  64 

16163 

14 

78.  1 

.  10 

310  Std 

1409 

28  May  64 

16164 

GT -60 

40 

.  50 

97.  1  Std 

1244 

28  May  64 

16171 

7 

78.  1 

.  50 

310  Std 

1)32 

1  June  64 

16172 

GT  - 14 

36 

.  10 

85.  1  Std 

144) 

1  J  une  64 

1617) 

46 

78.  1 

.  10 

310  Std 

1610 

1  June  64 

16174 

GT-11 

36 

.  10 

89.  9  Std 

1  327 

1  June  64 

Table  9 

Comparison  of  Variance*  for  GT  Balloon* 


Alt 

Alt 

Alt 

Alt 

GT  Balloons 

Std. 

0-10,000  11 

-20. 000 

20-30.000 

31 -40. 000 

16141  (GT  62-11) 

16144 

0.  037 

0.  942 

0.  894 

0.  739 

•16142  (GT  59-9) 

16144 

0.  812 

0.  947 

0.  954 

0.  842 

16146  (GT  56-6) 

16145 

0.  834 

0.  950 

0.  896 

0.  810 

16147  (GT  52-2) 

16145 

0.928 

0.  953 

0.  940 

0.  691 

16151  (GT  61-12) 

16152 

0.  892 

0.  903 

0.  817 

0.  740 

16154  (GT  58-8) 

161  57 

0.  890 

0.  965 

0.  961 

0.  91  5 

16155  (GT  55-5) 

16152 

0.  771 

0.  891 

0.  821 

0.  870 

16156  (Gl-  53-3) 

161  53 

0.  894 

0.  963 

0.  972 

0.  904 

16161  (GT  57-10) 

16163 

0.  733 

0.  945 

0.  958 

0.  956 

16164  (GT  60-7) 

16162 

0.  863 

0.  920 

0.  960 

0.  934 

16172  (GT  54-4) 

161  71 

0.  784 

0.  926 

0.  765 

N.  G. 

•  Indicate*  the  type 

that  I* 

tentatively  chosen 

a*  the  new 

standard  ROSE 

balloon. 


Table  10 

Comparison  of  Variance*  for  GT  Balloon*  Using  Linear  Wind* 


Alt 

Alt 

Alt 

Alt 

GT  Balloon* 

Std 

0-10, 000 

11-20. 000 

21-30, 000 

31  -40. 000 

16141  (GT  62-11) 

16144 

0.94  5 

0.  953 

0.  929 

0.  769 

16142  (GT  59-9) 

16144 

0.  770 

0.  929 

0.  965 

0.  879 

16146  (GT  56-6) 

16145 

0.  988 

0.  963 

0.  957 

0,  860 

16147  (GT  52-2) 

16145 

0.  940 

0.  96 2 

0.  950 

0.  784 
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7. 


GENERAL  CONCLUSIONS 


The  highest  degree  polynomial  necessary  to  describe  the  motion  of 
the  ROSE  balloon  is  in  almost  all  instances  the  first  degree.  If  the  data 
contains  extraneous  points,  they  should  either  be  eliminated  before  fitting 
a  linear  function  or  the  method  of  finite -diffe rences  can  be  used--the  latter 
being  faster  and  more  economical. 

The  Power  Spectral  Density  Analysis  Technique  when  applied  to  the 
balloon  data  produced  no  useable  results.  It  was  therefore  suggested  that 
physical  methods  of  reducing  the  oscillations  be  pursued  rather  than  trying 
to  remove  them  mathematically  from  the  radar  data. 

Some  of  the  modified  ROSE  balloons  show  significant  reductions  of 
self-induced  motion.  Extreme  care  must  be  used  in  selecting  the  type  of 
modification  to  be  applied  to  the  remaining  standard  ROSE  balloons.  With 
the  proper  modifications,  the  standard  ROSE  balloon  can  be  used  to  measure 
the  detail  winds  with  fair  reliability. 

The  smaller  lightweight  balloons  seem  to  reduce  wind  variance  most 
consistently.  They  also  have  better  wind  response  and  are  much  simpler 
to  manufacture  than  the  modified  ROSE  balloons. 


8.  FUTURE  WORK 

Future  work  will  be  in  two  areas:  The  first  will  be  a  more  complete 
analysis  of  the  modified  ROSE  balloons  selected  by  scientists  at  AFCRL  to 
be  flown  for  a  study  which  would  determine  a  typical  wind  profile.  The 
second  will  be  to  develop  a  more  convenient  data -handling  system  for  the 
new  lightweight  ROSE  balloons.  Such  a  system  will  permit  computations 
at  the  sounding  site  and  will  eliminate  the  delay  involved  in  data  reduction  at 
a  remote  location. 
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APPENDIX  1 


INDUCED  MASS 


The  values  presented  by  Leviton  for  the  wind  response  errors  in  th^  ^  ^ 
sph  erical  balloon  technique  for  measuring  winds  were  questioned  by  Reed1  ' 
The  basis  for  the  difference  between  values  computed  by  Leviton  and  Reed  was 
attributed  to  the  fact  that  Leviton  assumed  the  apparent  mass  term  to  be  zero. 


A  basic  reference  in  aerodynamics,  "Theory  of  Flight"  by  R.  von  Mises, 
page  573,  provides  the  basis  for  the  additional  term  in  the  expression  for  mass. 


"In  this  argument,  as  well  as  in  the  preceding  section,  it  has  been  assumed 
that  the  air  reactions  on  a  moving  body  depend  on  the  instantaneous  state  of 
velocity  only,  not  on  the  accelerations  (and  higher  derivatives).  It  is  obvious 
that  this  can  be  only  an  approximation  and  that  some  influences  on  the  accelera¬ 
tion  must  exist.  The  theory  of  irrotational  flow  of  a  perfect  fluid  gives  a  certain 
answer  to  this  question.  According  to  this  theory,  a  body  moving  in  a  fluid 
originally  at  rest  behaves  like  a  body  of  increased  inertia:  There  is  a  term  of 
apparent  mass  to  be  added  to  its  real  mass  w/g,  the  reaction  of  the  surround¬ 
ing  fluid  would  be  taken  intoaccount  when  the  sphere  is  assumed  to  have  the  mass 


W 


+ 


2 

3 


IT 


pa 


3 


and  to  move  under  the  influence  of  the  other  forces  (weight,  etc.)  alone.  This 
includes  (for  the  case  of  a  sphere)  the  result  expressed  in  D'Alembert's 
paradox  (Sec.  IX.  3)  that  no  reaction  exists  if  the  motion  is  uniform." 


The  volume  of  radius  a  is  given  by 


V 


3 

it  a 


Solving  for 


3 

a 


give  s 


(69) 


(70) 


Substituting  in  the  expression  for  mass  by  von  Mises  gives 


M 


P 


3  •  V  B 
Tit 


simplifying 


M 


1 

7  ? 


v 


B 


(71) 


(72) 


The  term  ^  P  *  Vg  18  equal  to  one  half  the  mass  of  the  air  displacedby 
the  sphere.  The  equations  of  motion  must  therefore  include  the  mass  of  the 
balloon  and  the  mass  of  the  air  that  is  accelerated  as  though  it  were  part  of 
the  body. 


65 


APPENDIX  2 


AN  ANALYSIS  OF  THE  OSCILLATORY  MOTION 
OBSERVED  IN  THE  RISE  OF  A  SPHERICAL  BALLOON 


In  the  course  of  this  work  it  has  become  apparent  that  certain  assump¬ 
tions  must  be  made  if  the  analysis  is  to  be  mathematically  tractable. 

From  wind-tunnel  pressure  measurements  it  is  well  known  that  the 
fluid  flow  at  the  front  of  a  spherical  object  is  essentially  incompressible  and 
mviscid  for  Mach  Numbers  appreciably  less  than  one.  However,  at  the 
rear  of  a  sphere  the  wake  is  turbulent  or  at  least  rotational  for  Re  greater 
than  100-200.  In  the  case  of  ROSE,  Re  =  10^  after  the  first  0.5  second  or 
so  of  travel.  Because  of  this  turbulence,  we  must  resort  to  approximate 
methods  since  not  even  the  simplest  turbulent  motion  has  yet  lent  itself  to 
a  strict  mathematical  analysis. 


However,  a  certain  amount  of  information  may  be  obtained  by  a  consider¬ 
ation  of  the  supposedly  potential  flo\*  at  the  front  of  the  balloon.  In  doing  this, 
of  course,  no  information  as  to  the  causes  of  any  oscillatory  motion  is  obtained 
since  these  causes  are  tied  up  (apparently)  in  the  unstable  wake. 


What  is  done  to  overcome  this  lack  of  information  is  to  accept  the 
oscillation  as  a  fact  and  to  apply  it  as  a  boundary  condition  to  the  equations 
governing  the  flow  at  the  front  of  the  balloon.  Due  to  a  lack  of  experimental 
data,  the  precise  form  of  the  balloon's  oscillation  is  not  known;  it  has,  how¬ 
ever,  been  assumed  that  this  oscillation  is  in  the  form  of  a  helical  spiral. 


We  take  axes  at  the  center  of  the  balloon  as  shown: 


The  x-axis  is  horizontal,  and  the 
y-axis  is  vertical.  The  z-axis  is 
normal  to  the  helical  axis. 


The  angular  velocity  of  the  balloon  is  denoted  by  Q,  the  vertical 
velocity  of  rise  by  W,  and  the  radius  of  the  spiral  by  H,  The  radius  of  the 
balloon  is  R.  Hence,  the  axes  are  translating  with  velocity  W  and  rotating 
with  angular  velocity  n.  The  velocity  of  the  air  relative  to  axes  fixed  in 
space  is  given  by  the  gradient  of  the  velocity  potential,  U.  The  velocity  of 
the  fixed  axes  relative  to  the  moving  axes  is  -  W  -  Hx  (T  +  cT)  where  7  is 
the  radius  vector  from  the  moving  axes  to  the  point  in  question.  Hence,  the 
velocity  of  the  fluid  relative  to  the  moving  axes  (V)  is  given  by 


V  =  VU-W-fix(r  +  d) 


(73) 
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Since  we  have  assumed  the  flow  to  be  inviscid,  the  boundary  condition 
on  the  surface  of  the  balloon  is 


Vr  =  0  (74) 

where  Vr  =  the  radial  component  of  velocity  relative  to  balloon,  and  on  the 
surface  r  =  R.  The  above  vectors  have  the  components: 


w  = 

(0,  W, 

0) 

(75a) 

d  = 

(0,  0, 

d) 

(75b) 

n  = 

(0,  n, 

0) 

(75c) 

"r  = 

(x,  y. 

z) 

(75d) 

n  x 

(3  +  r ) 

=  (n(7  +  d),  o, 

fix  j 

(75e) 

We  use  spherical  coordinate 

s  defined 

by 

x  = 

r  sin  9  cos 

y  =  r  sin 

0  sin  4,,  z  =  r  cos  0. 

Now,  Vr 

0  on 

r  =  R,  so 

V  = 
r 

:  0  =  V 

r  u  -  [V  +  n  x 

(?  +  ?)] 

r 

at  r 

=  R. 

(76) 

We  have  the  general  relationship: 


v  =  v  sin  0  cos  d>  +  v  sin  0  sin  <4>  +  v  cos  9  (77) 

r  x  y  z 


Hence  V  U  = 

r  9r 


r  =  R 


=  ft(z  +  d)  sin  0  cos  4>  +  W  sin  0  sin  $  -i?  x  cos  0 


or 


8U 

— -  =  ftd  sin  9  cos  <$>  +  W  sin  9  sin  $ 


r  =  R 


(78) 

(79) 


This  last  condition  along  with  the  fact  that  U-*0  as  r— teomakes  up  the 
boundary  conditions  on  U.  The  velocity  potential  (U)  must  satisfy  Laplace's 
equation: 


V2  U  =  0 


(80) 


is 


A  solution  to  this  equation  which  satisfies  the  above  boundary  conditions 

U  =  -iy-  (A  sin  0  cos  <$>  +  B  sin  9  sin  <*>).  (81) 

r 


It  follows  that 


8U 

9r 


=  R 


(A  sin  9  cos  <>  +  D  sin  0  sin  <*>) 


=  Qd  sin  0  cos  tp  +  W  sin  0  sin  <p 


(82) 
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The  above  equality  must  hold  for  all  0  and  4>,  hence: 


A  =  fid  and 


B  =  W 


or  A  =  - 


1  3 

i8 


and  B  =  -  i  R3W 


Thus,  the  final  expression  for  the  potential  becomes 
1  R3 

U  *  *  2  *7  (ftd  sin  0  cos  4>  +  W  sin  0  sin  4>) 


(83) 


Two  equations  connecting  the  four  unknowns;  ft,  d,  0,  and  <$>  may  be 
obtained  by  finding  the  coordinates  of  the  stagnation  point.  We  do  this  by 

finding  V  and  V  at  r  =  R  and  equating  them  to  zero. 

0  4> 


v0  =  V0U  -  [w  +  ftx  (7  +  d)] 


0 


(84) 


In  general,  cos  0  cos  <$>  +  cos  0  sin  4>  ~  vz  sin  0  (85) 

Hence,  V„  =  —  -  ft(r  cos  0  +  d)  cos  0  cos  4,  -  W  cos  0  sin  4>  -ftr  sin20cos<*) 

(86) 


r  =  R 


V  =  V 
4,  4, 


=  -  —  (ftd  cos  0  cos  4>  i  W  cos  0  sin  <$>)  -ft  R  cos  $  (87) 

U  -  [w  f  fix  (r  +  cT  )J  (88) 


in  general,  cos  <t>  *  vx  sin  4> 


therefore  V  = 


1 


au 


4>  r  sin  0  8$ 


- —  -  W  cos  4>  -*  ft(r  cos  0  +  d)  sin  4> 


(V  )  =4-  f  ftd  sin  4>  -  W  cos  4>  "I  +  ftR  cos  0  sin  4> 

»  r  =  R  2  1  J 

We  now  equate  expression  (87)  and  (91)  to  zero  and  obtain: 
3  cos  0  (ftd  cos  4>  +  W  sin  4>)  +  2ft  R  cos  <>  =  0 

3(ftRd  sin  4>  -  W  cos  4>)  -f  2ftR  cos  0  sin  $  =  0 


(89) 

(90) 

(91) 


(92) 

(93) 
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These  two  equations  contain  the  four  unknowns:  fl  (balloon  angular 
velocity),  d  (radius  of  spiral),  and  0,  <p  (coordinates  of  stagnation  point). 

As  is  well-known,  however,  two  equations  in  four  unknowns  are  a 
rather  useless  bit  of  symbolism  in  themselves.  With  this  in  mind,  we 
come  to  the  main  problem  of  this  paper:  How  to  find  something  which 
occurs  at  the  stagnation  point,  and  how  to  express  this  something  in  the 
form  of  two  or  mere  analytic  expressions. 

We  have  at  our  disposal  another  hydrodynamic  equation,  namely  Euler's 
equation  of  motion  or,  in  its  first  integral  form,  Bernoulli's  equation.  For  a 
system  of  rotating  and  translating  axes  such  as  those  employed  in  this 
problem,  Bernoulli's  equation  takes  the  form: 

P  =  J  P  H2D2  0  pgy  -  j  p  V2  +  f(t)  (94) 


where  p  =  pressure,  p 
cente  r  of  balloon,  D  = 
axis  of  spiral. 

There  are  good  physical  reasons  for  believing  that  the  pressure  is 
maximum  at  the  stagnation  point  (where  V  s  0).  This  may  be  rationalized 
by  considering  that  the  stagnation  point  is  the  most  physically  significant 
point  on  the  top-half  of  the  balloon.  Since  the  pressure  has  a  maximum  at 
some  point  on  the  balloon,  we  would  expect  that  this  maximum  would  occur 
at  the  stagnation  point.  Furthermore,  it  is  generally  known  that  in  the  case 
of  any  solid  body  moving  in  an  airstream  with  constant  velocity  the  stagna¬ 
tion  point  and  pressure  maximum  occur  at  the  same  point.  With  this  bit  of 
justification  behind  us,  we  now  proceed  to  maximize  the  pressure; 

The  pressure  is  a  maximum  if 


(,d  TT;  xJ) 


aoove 


)  =  tunction  oi  time,  y  =  height 

=  perpendicular  distance  from 


0  on  r 


R. 


Substituting  for  D  and  y  in  (3)  we  get 


p  =  y  pi}2(d2  +  2  rd  cos  0  +  r2  cos2  0  +  r2  sin20  cos2  <£) 


(95) 


-p  gr  sin  9  sin  4>  -  y  p  V2  +  f(t) 


dp 

w 


dp 


l  Z  2  Z  Z 

j  p  (]  (-2  rd  sin  0  -2 r^  sin  0  cos  0  +  2r  sin  0  cos  0  cos  $)  (9b) 


-  p  gr  cos  0  sin  <J> 
2  .  ,  2  2 


=  j  pn  (  -2r  sin  0  cos  $  sin  <*>)  -p  gr  sin  0  cos  6 


(97) 
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(98) 


These  last  two  expressions,  when  evaluated  at  r  =  R  and  equated  to  zero, 

g,VC  2  2 

U  sin  0  (R  cos  0  sin  4>  +  d)  +  g  cos  0  sin  4>  =  0 

S22R  sin  0  sin  <♦>  +  g  =  0,  provided  sin  0  cos  <$>  4  0 


(99) 


We  now  have  at  our  disposal  four  equations;  namely,  92,  93,  98,  and 
99  in  the  four  unknowns  Cl,  d,  0,  and  It  is  possible,  in  principle,  to  solve 
these  four  equations  for  t‘ie  unknowns. 


Equations  93,  98,  and  99  are  easily  solved  by  eliminating  0  and  <$>.  However, 
the  inevitable  conclusion  reached  is  that 

S2d2  =  0  (100) 


and  upon  further  analysis  we  see  that  this  implies  that  Cl  i  0,  but  that  d  =  0. 
This  is  interpreted  as  meaning  that  the  balloon  is  rotating  about  its  axis  and 
rising  vertically.  We  can  only  conclude  from  this  that  (1)  a  rising  balloon 
cannot  spiral  or  (2)  the  assumption  made  about  the  pressure  at  the  stagnation 
point  is  false.  If  we  accept  conclusion  (2)  as  more  likely  and  shelve  our  con¬ 
ceptual  misgivings,  we  can  search  for  another  set  of  relationships. 

If  the  pressure  is  not  a  maximum  at  the  stagnation  point  as  defined  above, 
where  is  it  a  maximum?  There  is  one  other  physically  significant  point  on 
the  upper  half  of  the  balloon.  This  point,  which  we  might  call  the  null- 
disturbance  point,  is  that  point  at  which  the  hydrodynamic  disturbance  in 
the  air  due  to  the  oalloon  is  zero.  That  is  to  say,  it  is  the  point  at  which  the 
air  motion  relative  to  the  balloon  axes  is  purely  helical.  The  fact  that  the 
disturbance  is  zero  at  this  point  is  expressed  mathematically  by  the  equation 


on  r 


vu  =  o  (1 01) 

This  implies  that  =  ^’U  =  VU  =  0,  and  that  V  U  *  0  identically 

=  R.  U  *  r  r 

V  U  =  ”  y  cos  0  (f2d  cos  <♦>  +  W  sm  <*>)  (102) 

V  U  =  7  (Od  sin  $  -  W  cos  $)  (103) 


The  above  two  expressions  when  equated  to  zero  yield; 

cos  0  ( f2d  cos  <*>  +  W  sin  4>)  =  0  and  (104a) 

f2d  sin  4>  -  W  cos  <$>  =  0  (104b) 

These  two  equations  will  be  supplemented  by  the  two  pressure  conditions 
used  before,  namely 


3p  _  9p 
30  ”  9ct> 


0 


at  the  null-disturbance  point. 
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The  pressure  is  given  as  before  by 


1  _2  .  ,2  ,  ,  Q  ,  2  2  A  t  2  .  2  2 
p  =  y  p  ft  (d  =  2dr  cos  9  4  r  cos  9  +  r  sin  9  cos  <$>) 


(105) 


-p  gr  sin  0  sin  $  -  j  p  +  f(t) 


^  =  pft^  (-dr  sin  9  -  r^  sin  9  cos  9  4-  rC  sin  9  cos  9  cos2 * 4  $)  (10b) 


-p  gr  cos  9  sin  4>  - 


QV 

W 


and 


3P  _ 


.2  2  .2.  . 


tt  .  aV 


•gE  =  -  p  ft  r  sin  9  sin  $  cos  <*>  -  p  f  -  sin  9  cos  <$>  -  p  V  * 


A  little  algrebra  gives; 

V-  •  ~  ~  \  ft^R^  sin  0  cos^(t»  -  2*  Wft  R  sin  9  *in  4>  co®  4> 


and 


2  2.  .2 
ft  R  sin  9  cos  9  sin  4> 


V  •  ^  WOR  cos  9  cos^  <♦>  -ft^R^  sin  4>  cos  4> 

3  .2  2  2  2. 

4  —  WAR  cos  0  sin  4>  4  ft  R  cos  9  sin  <j>  cos  <$> 


(107) 


(108) 


(109) 


Ihese  substituted  into  equations  106  and  107  give,  along  with  the  condition  that 

E  =  E  =  0,  that 

0  4> 


2  3  2  2 

-ft  d  sin  0  -  g  cos  0  sin  <$>  +  g  ft  d  sin  0  cos  <$> 

4  2  Wft  sin  0  8 in  <>  cos  $  =  0 


(110a) 


g  sin  9  cos  <$>  4  ^  Wft  cos  0  =  0 


(l  l  0b) 


The  preceding  four  equations  (104a,  104b,  110a,  and  110b)  are  analogous 
to  92,  93,  98,  and  99  of  the  previous  analysis  although  it  is  apparent  there  is 
little  similarity  between  the  forms  of  the  two  sets  of  equations. 

Liquations  104a  and  104b  yield  one  physically  significant  solution,  namely, 

W 

cos  0  =  0  and  tan  4>  =  r— ? 
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(Ill) 


These  two  conditions  when  substituted  in  equations  l  1  Oa  and  110b  yield 

«2d  =  0 

which  is  identical  with  the  result  obtained  under  the  assumption  that  the 
pressure  is  maximum  with  V  =  0. 

Thus,  two  essentially  different  sets  of  assumptions  lead  to  the  same 
physically  insignificant  result- -namely  that  the  balloon  does  not  spiral.  We 
are  again  left  with  the  two  alternate  conclusions  mentioned  previously: 

(1)  that  the  balloon  can't  spiral  and  (2)  that  the  pressure  assumption  is  false. 

If  we  again  accept  conclus ion  (2)  as  correct,  we  are  once  more  left  with 
but  two  valid  equations  (either  92  and  93  or  104a  and  104b),  and  the  problem 
again  is  to  find  two  more.  This  could  easily  be  done  if  an  expression  could 
be  found  for  the  pressure  in  the  vicinity  of  the  stagnation  point  (or  null- 
disturbance  point).  First,  we  recall  Bernoulli's  equation  is 

p  -  j  p  fl2D2  +  p gy  +  j  p  v2  =  f(t)  (ll^) 

The  above-mentioned  expression  for  the  pressure  could  be  substituted 
in  Bernoulli's  equation,  and  the  entire  system  of  equations  could  be  differ¬ 
entiated  with  respect  to  the  space  coordinates  of  the  balloon  and  equated  to 
zero. 


This  method  depends  upon  the  fact  that  f(t)  is  only  a  function  of  t  and 
the  refore: 


8f  _  8f  _  8f  92f  _  82f 

M  ~  9$  8r  90  8<t>  -  90  9r 


(113) 


Although  all  the  above  equations  are  not  independent,  at  least  three  of  them 
would  most  likely  be  so.  It  is  therefore  possible  that  they  could  yield  at 
least  three  more  equations  in  the  unknowns. 
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APPENDIX  3 


TIME  SERIES  ANALYSIS  BY  HEALY  AND  BOGERT 


1.  READIN  (A,  NA,  ANAME) 

This  Subroutine  reads  in  NA  values  of  the  sequence  A  from  punch  cards. 
NA  is  read  from  columns  1-6  of  the  first  card  and  the  BCD  characters  in 
columns  7-80  are  stored  in  the  vector  ANAME  for  future  references.  ANAME 
must  be  given  the  dimension  13  in  the  main  program.  The  second  card  of  the 
deck  carries  the  format  of  the  data  cards  in  ordinary  FORTRAN  form  omitting 
the  word  FORMAT  but  including  the  parentheses. 

2.  OUTPUT  (A,  NA,  ANAME,  INDIC) 

This  Subroutine  prints  NA  values  of  the  sequence  A  in  floating  form 
(without  starting  a  new  page).  The  printout  is  preceded  by  the  name  held  in 
BCD  characters  in  the  vector  ANAME  which  must  be  given  the  dimension  13 
in  the  main  program.  If  INDIC  is  negative,  the  series  will  be  punched  out 
five  to  a  card  with  the  two  preceding  cards  required  by  READIN;  if  INDIC 
is  a  positive  number,  the  series  will  be  written  on  the  corresponding  tape 
as  three  binary  records--the  first  containing  NA,  the  second  ANAME. 

3.  DETRND  (A,  NA,  B,  NDEG) 

The  mean  (NDEG  =  0)  or  a  least-square  linear  trend  (NDEG  =  1)  is 
subtracted  from  the  NA  values  of  the  sequence  A,  and  the  residuals  are 
stored  as  the  sequence  B. 

4.  AUTCOV  (A,  NA,  B,  L) 

This  Subroutine  calculates  the  series  B  as  the  autocovariances  of  the 
series  A  (of  length  NA)  from  lags  0  to  L.  Owing  to  the  nature  of  FORTRAN 
indexing,  B(J)  corresponds  to  lag  (J  -  1),  and  the  sequence  B  is  of  length 
L  +  1.  The  formula  used  is 

NA  -  J  +  l  2 

B(J)  =  L  A(I)  *  A(I  +  J  -  1)  /  (NA  -  J  +  1)  -  AVE  ,  (114) 

1=1  J  =  1,  L  +  1 

when  AVE  denotes  the  average  values  of  the  sequence  A. 

5.  FOUR  TR  (A,  L,  B,  INDIC) 

Form  the  sequence  B  as  the  Fourier  transform  of  the  sequence  A  of 
length  L  +  1.  For  INDIC  =  1, 

R(K)  =  A(1 )  4  £  A(I)  »  cos  (I  '  11  '  l)  it  4  (- 1  )K"  1  A(L+ 1 )  /  L,  (11S) 

1  =  2  K  =  l.  L  +1. 
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For  INDIC  =  2  the  cosines  are  replaced  by  sines,  and  the  end  terms  are 
omitted.  INDIC  =  3  or  4  provides  these  same  transforms  smoothed 
("hanned")  with  coefficients  1/4,  1/2,  1/4,  the  end  terms  being  found 
from  symmetry  considerations.  Note  that  the  sequences  are  not  of  length 
L+l  to  match  AUTCOV  and  CRSCOV.  N  must  not  exceed  1500. 
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APPENDIX  4 


ADDITIONAL  SUBROUTINES  ADDED  BY  WILLIAM  E.  BROCKMAN 
TO  THE  ORIGINAL  PROGRAM  BY  HEALY  AND  BOGERT 

1.  QTRND  (A,  NA,  B) 

The  least-squares  quadratic  trend  is  subtracted  from  the  NA  values  of 
the  sequence  A,  and  the  residuals  are  stored  as  the  sequence  B.  Orthogonal 
polynomial  coefficients  are  used. 

2.  WHITN  (A,  NA,  ALFAl ,  ALFA2,  B) 

The  data  in  the  sequence  A  of  length  NA  is  prewhitened  using  the  coeffi¬ 
cients  ALFAl  and  ALFA2;  the  whitened  data  is  then  stored  in  sequence  B. 

The  prewhitening  equation  is  given  by: 

B(I)  =  ALFAl  *  A(I)  +  ALFA2*A(I  +  1)  1  =  1,  N-l  (116) 

3.  PDRKN  (A,  NA,  ALFAl,  ALFA 2,  B) 

The  data  in  the  sequence  A  of  length  NA  is  post-darkened  (also  known 
as  "recoloring")  using  the  coefficients  ALFAl  and  ALFA2;  the  post -da rkened 
data  is  then  stored  in  sequence  B.  The  pos t -da rkening  equation  used  is 


B(I)  = 


A(I) 


ALFAl2  +  ALFA22  +  2*  ALFAl  *  ALFA2*cos  (tt*(I-1)  /  N) 

1=1,  N-l 


(117) 


4.  TRNDR  (A,  NA,  L,  B) 

Straight  line  segments  containing  L  points  are  subtracted  from  corre¬ 
sponding  points  in  the  sequence  A  of  length  NA;  the  residuals  are  stored  as  the 
sequence  B. 

5.  MOVAVG  (A,  NA,  L,  B) 

Moving  averages  over  L  points  are  subtracted  from  the  point  corre¬ 
sponding  to  the  mid-point  of  the  L  points  of  sequence  A  of  length  NA.  Again 
the  residuals  are  stored  as  the  sequence  B. 
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APPENDIX  5 


APPROXIMATE  VALUES  OF  AERODYNAMIC  FORCES 
ACTING  ON  THE  ROSE  BALLOON 


The  results  of  our  work  thus  far  in  explaining  mathematically  the 
balloon's  self-induced  motion  have  yielded  very  little  useful  information. 
The  plots  of  the  space  coordinates  of  a  ROSE  sounding  which  first  revealed 
loops  and  spirals  in  an  altitude  band  which  rawinsonde  reported  to  be  calm 
were  indeed  surprising.  Curiosity  was  the  motivation  for  our  attempt  to 
compute  the  magnitude  of  the  forces  necessary  to  produce  such  motions. 


The  equations  of  motion  for  the  three  coordinates  can  be  written  as: 

Mx  =  1/2  p  •  A  •  CD  |  v  |  (x  -  Wx)  +  Ax  (ll8) 

My  =  1/2  p  •  A  *  CD  |  v|  (y  -  Wy)  +  Ay  (119) 

M’i  =  -1/2  p  •  A  •  CD  I  v |  (z  -  W?)  +  A?  +  mg  -  gpVfi  ( 1 Z 0) 

M  is  the  total  (system  plus  induced)  mass,  and  A  is  the  Aerodynamic 
Force  in  the  q  component.  We  can  assume  that  Wz  =  z  =  A  =0,  and  we 
have  shown  that  a  good  approximation  for  M  (the  total  mass)  is 

M  =  m  +  1 /2  p  ( 1 Z 1 ) 


where  m  is  the  mass  of  the  system. 

Making  the  substitutions  and  rearranging  the  equations  become: 


-1/2  P  '  A  •  CD  |v| 


(m-i  1  /2  p  V 

■  ■  ■■  r 

X  - 


A 

x 


(122) 


-1/2  p  •  A  •  CD  |v|  = 

(m  +  1  /2p  Vfi)  y  -  Ay 

(123) 

y  -  wy 

-1/2  p  •  A  •  CD  |v|  = 

-mg  +  gp  Vfi 

(124) 

• 

7 

Since  the  left  side  of  all  three  equations  are  now  identical,  the  right  sides  can 
be  equated: 


(m+ 1/2  p  Vfi)  x  -  Ax 
x  -  W 

x 


-mg  +  gp  VB 
z 


( 1 2  5a ) 
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(125b) 


(m  +  1/2  p  V B)y  -  -  rng  +  gp 


and  solving  those  2  equations  for  A  and  A  respectively  gives: 

x  y 

(x  -  W  ) 

Ax  =  (m  +  1/2  p  Vfi)  x  +  (mg  -  gp  Vg)  •  - - — —  (12ba) 

(y  -  W  ) 

A  =  (m  +  l  /2  p  VQ)  y  +  (mg  -  gp  VQ)  •  - : — (126b) 

*  z 

The  values  used  for  the  accelerations  and  velocities  were  obtained  from 
a  randomly  selected  standard  ROSE  flight.  The  data  was  fitted  with  a  second- 
degree  orthogonal  polynomial.  The  first  derivative  evaluated  at  the  mid¬ 
point  of  the  100-foot  interval  produced  the  velocities,  and  the  second 
derivative  produced  the  accelerations.  The  density  (p)  was  obtained  from  the 
1962  Model  Atmosphere,  and  the  gravity  (g)  was  taken  to  be  -9.8  m/sec  . 

The  winds  (\V  x  and  Wy)  were  average  value  of  the  winds  on  the  layer  500  feet 
above  and  below  the  particular  altitude.  The  vector  sum  of  Ax  and  Ay  was 
then  computed.  The  results  were  as  follows: 

Altitude  (feet)  Aerodynamic  Force  (newtons) 


500 

27.  5 

5,  000 

45.  1 

10, 000 

28.  1 

20,  000 

10.  4 

30, 000 

11.5 

40, 000 

3.2 

50,  000 

5.5 

Ihe  values  seem  to  agree  with  the  spirals  in  the  data,  but  tneir 
accuracy  is  suspect  because  of  approximations  and  assumptions.  The  trend 
of  lower  forces  at  the  higher  altitudes  also  agrees  with  the  theory  that  the 
critical  value  of  Reynolds  Number  occurs  at  about  40,000  feet. 

This  approach  was  not  pursued  because  of  reasoning  similar  to  that 
given  for  not  pursuing  the  questionable  results  of  Power  Spectral  Density 
Analysis. 
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APPENDIX  6 


LEGENDRE  POLYNOMIALS 


This  appendix  is  included  in  order  to  give  the  reader  easy  access  to  the 
aitu.il  numbers  used  in  the  polynomials.  Pm  n(i)  is  a  Legendre  polynomial 
of  degree  m,  of  the  form 


m  .  .  [k] 

pr  n(i)  -  E  (VO  1 

m,n  k  =  0  k  k  [k] 

n  J 

(13) 

w  he  re  i  ^k^  =  i  (i -  I  )(i -2  )  .  .  .  (i  - k-t  1 )  and  i  =  1 

(14a) 

and 

.  m .  ml 

(k'  "  TTfTkTTTT 

(14b) 

The  first  few  polynomials  are; 

P0,n<'>  =  ‘ 

(33a) 

Pl.n<‘>  =  1  -  ^ 

(33b) 

p  (i)  =  1  -  6-  3  6  'j"!) 

2,  n  n  n(n- 1 ) 

(33c) 

Now 

make  the  following  substitution; 

C1  =  - 
1  n 

(127a) 

C2  *  C1  <  1  -T  > 

(127b) 

c3  *  S  <  hr » 

• 

(127c) 

• 

C  =  c  1  (  — 7T  ) 
m  m-ln-m-tl 

(1  27d) 

give  s 

this  simpler  form  for  the  polynomie  Is . 

P0.n<»  *  1 

(128a) 

Pl,n<‘»  *  1  -  iCl 

(128b) 

P2.n(i>  *  1  -  6C1  *  6C2 

(128c  ) 
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P?n  (i)  =  1  -  IZC{  +  30  C2  -  20C3  (128d) 

P4  n(i)  =  1  -  20CJ  +  90  C2  -  140C3  +  70C4  (128e) 

Pc  (1)  =  1  -  30C,  +  210C,  -  560C,  +  630C  .  -  252Cc  (128f) 

5,  n'  '  1  2  3  4  5 

Pt  (1)  =  1  -  42C.  +  420C,  -  1680C,  +  3150C„  -  2772Cc  +  924C,  (128g) 

0,  n  I  L  3  4  3  0 


P,  (1)  =  I  -  56C,  +  756C,  -  4200C,  +  1  1550C.  -  16632C.+  12012C,-  3432C-, 
7-n  1  2  3  4  5  6  (128h)7 

PQ  (i)  =  1-  72C.  +  1260C,  -  9240C, +34650C,  -  72072Cc  +  64084C, -51480C., 

8,  n  1  2  3  4  5  b  7 

+  1 2870Cg  (12  81) 

Prt  (1)  =  1  -  90C,  +  1 980C,  -  18480C,  +  90090r  -  2 52252Cc  +  420420C, 

9, n  1  2  3  4  5  6 

-  411840C.  +  218790Co  -  48620C_  (1Z8J) 

7  07 


The  first  derivative  of  the  orthogonal  polynomial 


■ 


A0.n<‘> 


AlPl.n<‘» 


♦A2P2.n“> 


A  P  (i) 
m  m,  n' 


(11) 


is  given  by 


-  Aop,o.n<'> 


Aip,l.n<l> 


A2p,2.n(i> 


+  A  P*  (i)  (30) 

m  m,  n 


\«.here  P1^  ^(i)  is  the  first  derivative  of  the  Legendre  polynomial 

It  is  a  characteristic  of  Legendre  polynomials  that  if  m  is  even  and  the 
first  derivative  is  evaluated  at  the  mid-point  (n/2)  of  the  interval  with  length  n, 
then 


P 


n 


m  n{J)  =  ° 
m,  n  c 


(38) 


Therefore,  the  derivative  of  an  even  degree  polynomial  is  equivalent  to 
the  derivative  of  a  polynomial  of  one  degree  less.  For  example,  if  m  =  2, 
then 


(i>  =  1 


n 


6 


2 

1  -  1 

"Z - 

n  -  n 


(33c) 
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(34c) 


and  the  first  derivative  with  respect  to  i  is  given  by: 

P'  (i)  =  -  6  •  i  +  6  • 

2,  n  '  n  Z 

n  -  n 


At  i  =  n/Z 


6  +  6(Z  •  ?  -  n 
n  n  (n  -  1 ) 


( 1  Z  9a ) 


lb  b  (n-  1 
n  n  n  -  1 


6  6 
n  +  n 


(1  29b) 


The  derivatives  are  taken  with  respect  to  the  index  (i);  therefore  the 
denominator  which  contains  only  functions  of  n  does  not  change.  To  simplify, 
the  following  substitution  will  be  made:  Let 


NB  (1)  =  n 
NB  (Z)  =  n  (n - 1  ) 

NB  (3)  =  n  (n - 1 )  (n-Z) 


(130a) 
(130b) 
(1  30c) 


NB  (j)  =  n(n-l)  .  .  .  (n+  1  -  j) 


( 1  3  Od ) 


The  Legendre  polynomial  general  derivatives  and  the  derivatives  evaluated  at 
i  =  n/Z  are  listed  below: 


P'l,n(l)  =  -N5PT 


P'l,n(I)  =  ‘  5TFfTT 


p,3,n<'»  *  •  12  mw  3  30  K) 


(3i2  -  6i  +  Z) 


(131a) 

(131b) 


(132a) 


p,3,n(Z»^m  -  30N5iI,['’-1]-Z0NBT3)[2+Z«-,Zt3n)]  <132b> 


P'5,n“»=  -3°NF7T)+  21°N^TZt[2'  -  >]  *  560n1T3)[2  +  I  <*»2  +  <>>)]<133.l> 


+  630  nF(T)['  6  +  ,<22  3  2,(’9  +  2l,)]  ‘  252 


NB  (5)  * 


Z4  +  ^  (  -  800  4  Zi  (420  4  2i(  -  80  +  1  Oi)))J 


P'5,n<I>  '  30NBTU 


3  210  NF|T)[n  ‘ ']  '96°  n573)[2  3  ?  «'12  f  3">]<l33,'> 
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P' 


+  630 


i  r  n 

NB«)  L6  +  2(22  +n(-9  +  n|)]  -  252  1 

J  TO75)  ' 


L24  +  16  (  -  800  t  n(420  +  n(  -80  t  5n)))] 

pt  »n.  I 

7.  n  F  =  ~56  mr/,7  +  756  — 1  T  ,  ,  1  i 

NB(1)  NB  (2)  I  _1  +  n  -  4200  1 

1  '  L  J  NBlT)*  H34) 

[2  +  4  (-12  +  3n)l  +1150  _ L__  f  ,  n 

J  NB  (4)  L  “  6  +  F(2Z  +  n(“9  +  „))1 

-  16632  - L_  f p  ^  ,  n  J 

NB  (5)  L  ~M°0  +  "(420  +  n(_80  +  5n)))J 

+  12012  L_  T  x  n  /. 

NB(6,  l-lM+1-,4348to(.„00  +  nN75+jn 

-3432  —  ^  f  7,.  n 

NB(7,l  720  +64  <-*'2«96  +  n(77952  +  7, . 12000  +  n 
(3500  +  n(252  +  7n)))))J 

NB  +  1980  ffiffT,  f'1  +  "]  -  18480  _L_.  (135) 

f' *7'-“ •’"']* -R,  n] 

-252252  - - ! -  f  j.  n  / 

NB  (5)  *•  16  1_8°8  +  n{420  +  n{-80  +  5n),)"* 

+  420420  — ! -  I"  n  , 

NB,6,  +  “  (4384  +  n( -2700  +  n(-75  +  3n ) ) ))"* 

'41184°  NbTtT  [  720  +-£(-U2«96+n(77952  +n 

('12000  +  n(3M°  •  "<252  4  ?"»»>]  +  218790  _J__ 
r  NB  (8)  • 

[-8040  ♦  -  (209088  .  57584  ♦  n(54152  4  n(-980  +  „ 

-,.9  4  ■  »»„]  .48820  [40320  +-^L 


9  .n(F>  =-9° 


f~28053504  +  n(2267,808  +  a(-8612352  +  n(17,5,20  +  n- 
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(-217728  +  n(  1  5288  +  n(-576  +  9)))))))] 


A  still  more  convenient  form  for  computing  can  be  obtained  by  the  following 
substitution: 


Let  (FCR^  = 
(FCR)Z  = 

(fcr)3  = 

(FCR)4  = 
(FCR)_  = 


1 


1  NB  (1) 
1 

2  'NB  (2) 
1 


(FCRL  = 


3  NB  (3) 
1 

4  NB  (4) 

1 _ 

5  NB  (5) 

1 


-1  + 

-12  +  3 


■] 


■ 

n 


n 


■6  +  —  (22  +  n(  -9  +  n)) 

b 


n 


24  +  -£  (-800  +  n(420  +  n(-80  +  5n))) 
1  6 


(fcr)„  = 


6  NB  (6) 
1 


n 


-120  +—  (4384  +  n(-2700  +  n(680  -  n(-75  +  3n)))) 
1 6 


(136a) 

(136b) 

(136c) 

(1 36d) 

(1 36e) 

(13  6f ) 


7  NB  (7) 


n 


720  +  —  (-112896  +  n(77952  +  n(-12000  +  n  (1 36g) 

64 


(3500  +  n(252  +  7n))))) 


(FCR)0  = 


1 


8  NB  (8) 


n 


-5040  +  —  (209088  =  n(-157584  +  n(54152  +  n(-9800  (1  36h) 


+  n(966  +  n(-49  +  n)))))) 


(FCR)ft  = 


1 


9  NB  (9) 


n 


40320  +—T7  (-28054504  +  n(22679808  +  n(-8612352  (1  36i) 

COD 


+  n(  1 795920  +  n(-21  7728  +  n(  15288  +  n(-576  +  9n))))))) 

The  polynomials  evaluated  at  the  midpoint  (n/2)  then  become: 

P'.  £)  =  -2  (FCR).  (137a) 

1 ,  n  c  l 

P'  (t)  =  -12  (FCR),  +  30  (FCR),  -  20  (FCR),  (137b) 

3,  n  2  1  2  3 

P'  (7)  =  -30  (FCR)  +  210  (FCR)  -  560  (FCR)  +  630  (FCR)  (137c) 

D|  n  4  1  b  j  *t 


-252  (FCR), 
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p,7  (7> 

7,  n  2 


P 


9.  „<7> 


-56  (FCR),  +  756  (FCR)  -  4200  (FCR)  +  11550  (FCR)  ( 1  3 7d) 

12  3  4 

-16632  (FCR)  +  12012  (FCR),  -  3432  (FCR) 

5  6  7 

-90  (FCR),  +  1980  (FCR)  -  18480  (FCR),  =  90090  (FCR),  (137o) 
12  3  4 

-252252  (FCR),.  +  420420  (FCR)6  -  41  1840  (FCR)? 

+  218790  (FCR)  -  48620  (FCR) 

8  9 
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